arXiv:hep-ph/9502254vl 8 Feb 1995 


Centre de Physique Theorique*- CNRS - Luminy, Case 907 
F-13288 Marseille Cedex 9 - France 


CHIRAL LAGRANGIANS 

AND 

KAON CP-VIOLATION * 1 

Eduardo de RAFAEL 


Abstract 

These lectures are an introduction to the subject of chiral effective Lagrangians of 
the Standard Model and their applications, mostly in the sector of non-leptonic kaon 
decays, with special emphasis on CP-violation. The first lecture gives an introduction 
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1 Phenomenology of K°—K° Mixing and CP—Violation 

in K —> 2n Decays 

The purpose of this first lecture is to introduce the basics of the phenomenological 
description of the K°-K° system. There is practically no theory behind this de¬ 
scription. It is only based on first principles: the superposition principle, Lorentz 
invariance, and general invariance properties under the P, C and T symmetries. The 
basic idea is to reduce the description of this system to a minimum of phenomenolog¬ 
ical parameters which, eventually, an underlying theory -like the Standard Model- 
should be able to predict. There are many reviews on this subject. I recommend 
the reader to consult Ref. 0 for an excellent historical overview and complementary 
information. 

1.1 Phenomenology of K°—K° Mixing 

In the absence of the weak interactions, the K° and K° particles produced by the 
strong interactions are stable eigenstates of strangeness with eigenvalues ±1. In 
the presence of the weak interaction they become unstable. The states with an 
exponential time dependence law (r is the proper time) 

| K l ) - e~ iMLT | K l ) and | K s ) - e~ iMsr \ K s ), (1) 

are linear superpositions of the eigenstates of strangeness: 

\K l ) = 1 (p K°) + q I K»)) 

vl P I + I 9 I 

W = /, W=, Tf (p I K°) - V T°>) , 

VI p I +191 

where p and q are complex numbers and CPT-invariance, which is a property of the 
Standard Model in any case, has been assumed. The parameters M Lt s in Eq. (|T|) 
are also complex 

M LjS = m L ' S - 2 -T L ,s , ( 4 ) 

with rriL,s the masses and Tl^s the decay widths of the long-lived and short-lived 
neutral kaon states. 

As we shall see, experimentally, the | K$) and | Kj) states are very close to the 
CP-eigenstates 

with 

CP | Kl) = + | K°) and CP | A 2 °) = - | K° 2 ). (6) 


)+ I K»)) (5) 


K °} = ^{\ K 0 ) -\K a }) and \K°) 


( 2 ) 

( 3 ) 
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This is characterized by the small complex parameter e, 

~ P~q 
e = -; 

p + q 


( 7 ) 


in terms of which, 


Kps) 



(I K\ 


2 , 1 ) + 


K° h2 )) ■ 


( 8 ) 


According to Eqs. i) and (H), a state initially pure [ K°) evolves, in a period of 
time r to a state which is a superposition of | K°) and | K°) : 


| K°) -> + e ~ iMsT ] | K°) + \- q \ e ~ iMLT - e~ iMsT } \ K°) ; (9) 

and, likewise 

I K°) -»■ + e -* Msr ] | K°) + - e" iMsr ] | K°). (10) 

z z p 

For a small period of time dr we then have 


| K°) - 

> | K°) - iSr(M n 

| K°) +M 12 

K 0 )); 

(ii) 

1 K°) - 

> | K°) - iSr(M 21 

| K°) + M 22 

1 77°)), 

(12) 


where 




If M l + M s ^M l -M s )\ 

2 [ q -(M L -M s ) M l + M s )■ 


(13) 


This is the complex mass-matrix of the K° — K° system. 

In full generality, the mass-matrix Adp admits a decomposition, similar to the 
one of the complex parameters Ml,s in Eq. ([|). in terms of an absorptive part Tp 
and a dispersive part Mp : 

Mij = Mp - ^Tp . (14) 


In a given quantum field theory, like e.g., the Standard Electroweak Model, the com¬ 
plex K° — K° mass-matrix is defined via the transition matrix T which characterizes 
S'-matrix elements. More precisely, the off-diagonal absorptive matrix element T \ 2 
for example, is given by the sum of products of on-shell matrix elements: 


T 12 = J2 rfr (( r I T I K 0 ))* (r | T | K°) , (15) 

r 7 

where the sum is extended to all possible states | T) to which the states | K°) and 
| K°) can decay. The symbol dr denotes the phase space measure appropriate to 


2 





the particle content of the state T. The corresponding matrix element M \2 is defined 
by the dispersive principal part integral 

1 r 1 

M 12 = — p / ds — 5 -IWs) + “local — terms” . (16) 

7 r J rn K — s 

The fact that M.\\ = M .22 in Eq. (|T^) is a consequence of CPT-invariance. In 

general, if we have a transition between an initial state | IN) and a final state | FN ). 

CPT-invariance relates the matrix elements of this transition to the one between 

_ / _ / _ / 

the corresponding CPT-transformed states | FN ) and | IN ), where | IN ) denotes 
the state obtained from | IN) by interchanging all particles into antiparticles (this is 
the meaning of the bar symbol in IN), and taking the mirror image of the kinematic 
variables: [(E,p) —> (E, — p ); (cr°, a) —> (—a 0 , a)], as well as their motion reversal 
image: [(E,p) —> (E, — p) ; (cr°, a) —> (cr 0 ,—a)]. (These kinematic changes are the 
meaning of the prime symbol in IN ).) Altogether, CPT-invariance implies then: 

(FN \T | IN) = (In' \T I FN). (17) 

Since, for the Ji°-states: j (K 0 ) ) —\ 77°), the CPT-invariance relation implies 


M 11 — M.22 ■ 


(18) 


The off-diagonal matrix elements in Eq. (|I3|) are also related by CPT-invariance, 
plus the hermiticity property of the T-matrix in the absence of strong final state 
interactions; certainly the case when the | IN) and | FN) states are | K°) and 
| A' 0 ). In general, in the absence of strong final state interactions, we have: 


(IN | T | FN) = ({FN | T | IN)) 


(19) 


This relation, together with the CPT-invariance relation in (|17) implies then: 


M.12 — (M.21Y ■ 


( 20 ) 


There are a number of interesting constraints between the various phenomeno¬ 
logical parameters we have introduced. With M 12 and r 12 defined in Eqs. ([16]) and 
O and using Eqs. (]T^, [| and |7j), we have 


q 1 — e 1 Am + i|AT 

1 

to 

1 

tO|t^. 

to 

(21) 

p 1 + 7 2 M 12 — §ri2 

" 5 (Am + |Ar) ’ 

where 



Am = mi — ms and 

Ar = r s -r L . 

(22) 

As already discussed, CPT-invariance implies 



M 2 i = (A/ 12 )* and 

r 2 i = (ri 2 )*. 

(23) 
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Experimentally, the masses rriL,s and widths r l,s are well measured, and in 
what follows they will be used as known parameters. (There is no way for theory at 
present to do better than experiments in the determination of these parameters...) 
The precise values for the masses and widths can be found in the Particle Data 
Booklet a (PDB). Nevertheless, it is important to keep in mind some orders of 
magnitude: 


p-l 

1 s 

~ 0.9 x 10 10 sec.; 

(24) 

r L 

~ 1.7 x 10~ 3 r 5 ; 

(25) 

Am 

~ o. 5 r 5 . 

(26) 


1.2 The Bell—Steinberger Unitarity Constraint 

Let us consider a state | T) to be an arbitrary superposition of the short-lived and 
long-lived kaon states: 

|T) = a\K s )+P\K L ). (27) 

The total decay rate of this state must be compensated by a decrease of its norm: 


£ I (r I T I T) | 2 = - — I T I 2 . 

P C LI 


(28) 


The change in rate is governed by the mass matrix defined by Eq. (0)- Equating 
terms proportional to | a | 2 and | (3 j 2 in both sides of Eq. (^) results in the trivial 
relations: 


r L = J2 [ d r 1 (r 1 t | k l ) i 2 , 

r J 

(29) 

Y s = J2 f d r 1 (r 1T 1 Ks) I 2 • 

r ^ 

(30) 


The mixed terms, proportional to a/3* and a*(3. lead however to a highly non-trivial 
relation, first derived by Bell and Steinberger B: 

- KM* L - M S )(K L I Ks) = £ / dr «r m k l ))* (r | t \ k s ) . (31) 

r J 


Notice that 


{K l I K s ) 


1 P 

2 — 

q 

| 2 2Ree 

1 P 1 

2 + 

q 

| 2 1+ 

e 

2 


(32) 


The l.li.s. of Eq. ( fTTl ) can be expressed in terms of measurable physical parameters 
with the result 


(rj+qi __|R«_ = Y^jdv «r | t | k l ))' (r | t k s ). (33) 
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The r.h.s. of this equation can be bounded, using the Schwartz inequality, with the 
result 


r 5 + r L 

2 


— iAm 


2 Ree 

1+ I e I 2 


< i/TlTs ■ 


(34) 


Inserting the experimental values for T s,l and Am, results in an interesting bound 
for the non-orthogonality of the Kl and Ks states [see eq ([32P1: 


2 Ree 
1 + I e I 2 


< 2.9 x 10~ 2 , 


(35) 


indicating also that the admixture of K®(K%) in Kl(K s) has to be rather small. 


It is possible to obtain further information from the unitarity constraint in ([33]), 
if one uses the experimental fact that the 2tt states are by far the dominant terms 
in the sum over hadronic states T. One can then write the r.h.s. of Eq. ([33]) in the 
form 

E f h(vr 7 r) ((7T7T I T I K l ))* (tttt I T I K s ) + yr 5 . (36) 


It is possible to obtain a bound for 7 , by considering other states than 2n in the sum 
of the r.h.s. in Eq. ([33]) and applying the Schwartz inequality to individual sets of 
states separated by selection rules. The contribution from the various semileptonic 
modes, for example, is known to be smaller than 


| £ /••• |«10- 3 r 5 ; (37) 

lep. modes 

and the contribution from the 37 T-states 

IE /•••l«io- 3 r 5 . (38) 

3tt J 

We conclude that to a good approximation we can restrict the Bell-Steinberger 
relation to 27r-states. We shall later come back to this inequality, but first we have 
to discuss the phenomenology of the dominant K —> tttt transitions. 


1.3 K — > tttt Amplitudes 

In the limit where CP is conserved the states Ks(Kl) become eigenstates of CP; i.e., 
the states K^K®) introduced in (|5|) with eigenvalues CP = +1(CP = — 1). O 11 the 
other hand a state of two-pions with total angular momentum J = 0 has CP = +1. 
Therefore, the observation of a transition from the long-lived component of the 
neutral kaon system to a two-pion final state is evidence for CP-violation. The 
first observation of such a transition to the 7 r + 7 r - mode was made by Christenson, 
Cronin, Fitch, and Turlay □ in 1964, with the result 


r A+,~) 

r L (all) 


(2 ±0.4) x 10" 3 . 


(39) 
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Since then the transition to the 7r°7r° mode has also been observed, as well as the 
phases of the amplitude ratios 


V+- 


(7T + 7T 

T 

\K l ) 

(7T+7T - 

T 

Ks) 


with the results 


I: 


and 


Voo 


(7r°7T 0 I T I K l ) 
(7T 0 7r 0 j T I K s ) ’ 


?/+_ = (2.269 ± 0.023) x 10 _3 e i(44 - 3±a8) ° ; 
7700 = (2.259 ± 0.023) x 10“ 3 e i(43 ' 3±1 ' 3) ° . 


(40) 


(41) 

(42) 


In order to make a phenomenological analysis of K —> 7T7t transitions, it is 
convenient to express the states | 7r + 7r _ ) and | 7r°7r°) in terms of well dehned isospin 
I = 0, and 1 = 2 states. (The 1 = 1 state in this case is forbidden by Bose 
statistics.): 

+-} = \j^ I °> + ]f^ I 2); (43) 

00> = vf |2> -\/| |0> ' (44) 

The reason for introducing pure isospin states, is that the matrix elements of tran¬ 
sitions from K° and the K° states to the same (7T7 t)/- state can be related by CPT - 
invariance plus Watson’s theorem on final state interactions. The relation in ques¬ 
tion is the following 


e~ 2iSl (I | T | K°) = ((/ | T | K°))* , (45) 

where 5/ denotes the appropriate J = 0, isospin / tttt phase-shift at the energy of 
the neutral kaon mass. 

The proof of this relation is rather simple. With 5 = 1 + iT , the unitarity of 
the 5-matrix, 55^ = 1, implies 


T t T = i(T t -T). (46) 

If one takes matrix elements of this operator relation between an initial state K°, 
and a final 2x -stat e with isospin /, we then have 

J2(I \T ] \F)(F\T\ K°) = i(I | T f | K°) — i(I \ T \ K°), (47) 

F 

where we have inserted a complete set of states J2 \ F)(F \= 1 between T and Tk 
The crucial observation is that, in the strong interaction sector of the 5-matrix, 
only the state F = I can contribute to the Y' mat rix element. All the other states 
are suppressed by selection rules; e.g., the 37r-states have opposite G parity than 
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the 27r-states; the nlv- states are not related to 27r-states by the strong interactions 
alone; etc. Then, introducing the mr phase-shift definition: 

(/ | S | /) = e 2iSl , (48) 

results in the relation 

i(e~ 2iSl -1)(I\T\ K°) = i(I | T f | K°) — i(I \ T \ K°), 

= i((K° | T | I))* — i(I | T | K °). (49) 

We can next use CPT-invariance [recall Eq. (]T7[) . which in our case implies the 
relation: (. K° \ T \ I)* — ((/ | T \ Ji 0 ))*.] The result in Eq. (|45| ) then follows. 

As a consequence of the relation we have proved, we can use in full generality 

the following parametrization for K°(K°) (nn): amplitudes: 

(I \T\ K°) = iA ie ibl ; (50) 

(/ j T | K°) = -iA*e i51 . (51) 


One possible quantity we can introduce to characterize the amount of CP- 
violation in K —> 27T transitions is the parameter 


£ _ A[K l -> (7T7r) J=0 ] 
A[K S -> (vr7r)/ =0 ] ' 


This parameter is related to the e-parameter introduced in Eq. (^)); as well as to 
the complex A 0 -amplitude defined in ( |T)| ) and (|5l]) , in the following way 


(1 + e)Ap — (1 — e)A* Q 
(1 + e)A 0 + (1 — e)Ag 


i.e., 


e + i 


e = 


ImAp 

ReAo 


1 + ie 


ImAp 

Re^o 


(54) 


This is a good place to comment on the history of phase conventions in neutral 
K decays. In their pioneering paper on the phenomenology of the K — K system, 
Wu and Yang Q chose to freeze the arbitrary relative phase between the K° and 
K° states, with the choice ImAo = 0. With this convention, e — e. In fact, the 
parameter e is phase-convention independent; while neither e, nor A/ are. Indeed, 
under a small arbitrary phase change of the K {] st at e: 


1 K°) - 

-> e~ ilf | K °), 

(55) 

the parameters A/, M 12 , and e change 

as follows: 


Irn A r —> 

ImA/ — (pReA/; 

(56) 

ImM 12 — > 

ImM 12 + tpAm; 

(57) 

e — > 

e + iip] 

(58) 
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while e remains invariant. The Wu-Yang phase convention was made prior to the 
development of the electroweak theory. In the standard model, the conventional 
way by which the freedom in the choice of relative phases of the quark-fields has 
been frozen, is not compatible with the Wu-Yang convention. Since e is conven¬ 
tion independent, we shall keep it as one of the fundamental parameters. Then, 
however, we need a second parameter which characterizes the amount of intrinsic 
CP-violat.ion specific to the K —> 2tt decay, by contrast to the CP-violation in the 
K° — K° mass-matrix. The parameter we are looking for has to be sensitive then 
to the lack of relative reality of the the two isospin amplitudes Ao and A 2 . This is 
the origin of the famous e -parameter, which we shall next discuss. 

In general, we can define three independent ratios of the K L S —> (2 tt) /=(k2 tran¬ 
sition amplitudes. One is the e-parameter in (|52|). Two other natural ratios are 


A[K L -> (7T7t) J=2 ] = A[K S -> (7T7t) J=2 ] 

A[K S -> (vT7r) /=0 ] A[K S -»■ (7T7 t) /=0 ] ' 


(59) 


Both ratios can be expressed in terms of the e-parameter introduced in Eq. and 
the complex ^/-amplitudes defined in ([50]) and (|5l|): 


and 


A[K l -»• (7T7r)/ =2 ] 

A[K S -> (7T7r)/ =0 ] 


(1 + e)A 2 - (1 - e)A^ c i(s 2 -s 0 ) 
(1 + c)Aq + (1 — c)Aq 


: ImA? 
' ReA 0 


+ e 


1 + ie 


ReAa 

ReAo i(52—<5o) . 


ImAn 

ReAo 


A[K S (7T7T)j =2 ] 
A[K S -> (7T7r)/ =0 ] 


(1 + e)A 2 + (1 - e)A* ^ i ( S? _5 ()) 
(1 + e)A 0 + (1 — c)Aq 

ReA 2 1 ~ ImAg 

ReAp _ RcAp ^,( 62 — 60 ) 

1 + 


(60) 


(61) 


The e'-parameter is then dehned as the following combination of these ratios : 


/ 

e 


1 

71 


/ A[K L —> (7T7T)J =2 ] 
\A[K S (7T7r)/ =0 ] 


e x to 


(62) 


From these results, and using the expression for e we obtained in Eq. (|Y3|), we finally 
get 

i (T — 

'' = 75 (Re^ + fe-ImA,)^ (lm/l2Re ' 4 ° " Im - 4 » Re ^) • < 63 ) 

an expression which clearly shows the proportionality to the lack of relative reality 
between the A 0 and A 2 amplitudes. 
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We shall next establish contact with the parameters rj + - and r/oo, which were 
introduced in Eq. ([40], and which are directly accessible to experiment. Using 
Eqs. (f|3]), (pH|), as well as the definitions of e, e , and u above, one finds 


V+- = e + e 


1 + V 2 UJ 


Voo = 


2 e 


1 - ypk 


w 


(64) 

(65) 


So far, we have made no approximations in our phenomenological analysis of 
the K° — K° mass-matrix and K —> 2tt decays. It is however useful to try to thin 
down in some way the exact expressions we have derived, by taking into account 
the relative size of the various phenomenological parameters which appear in the 
expressions above. The strategy will be to neglect first, terms which are products 
of CP-violation parameters. For example, in Eq. (|6lj), we have introduced the 
parameter u, which a priori we can reasonably expect to be dominated by the term 


u ~ ———— ^o). 

ReA n 


( 66 ) 


We can justify this approximation by the fact that non-leptonic A/ = | transitions, 
although suppressed with respect to the A/ = | transitions, are nevertheless larger 
than the observed CP-violation effects. Notice that the amplitude A 2 is responsible 
for the deviation from an exact A/ = I rule. The ratio can be obtained 

from the experimentally known branching ratios P (Kg —*■ 7r + 7r - ) and P ( Ks —> 
7 r°7r°). More precisely, correcting for the phase-space effects, one must compare the 
normalized decay rates: 


7(1.2) = 


r(/i - 7Ti7r 2 ) 


1 


I67 tM 


(m-i+m-2) 2 

M 2 


1 - 


(mi—m2) 2 

M 2 


1 - 


(67) 


where the denominator here is the two-body phase space factor for the mode K — 
7r 1 7r 2 , (M is the mass of the /P-particle and m 12 the pion masses.) Then, we have 


i§oo| “ 1 + 3 v/I s ^ cos( ' 52 " Sl)+ °fy‘ 


Experimentally, from the PDB @, one finds 

7s(+ 


27 s (00) 


= 1.109 ±0.012, 


( 68 ) 


(69) 


and using the present experimental information on {82 — 5 1 ) - [see the discussion in 
Sec. 5.,] we find, with neglect of radiative corrections 


RePE 

ReAr 


= (± 22 . 2 )” 1 . 


(70) 
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I shall later discuss some of the qualitative dynamical explanations, within the 
standard model, of how this small number appears. It is fair to say however, that 
a reliable calculation of this ratio is still lacking at present. 

Using the approximations 


e ImA 0 Rev4 0 

we can rewrite e' in a simpler form 


and 


e 2 < 1 , 


ImA 


e ' ~ o+5) ( l ^l _: 

y/2 ReA 0 v ReA 2 ReAo 


clearly showing the fact that e' is proportional to direct CP-violation in K 
transitions and is also suppressed by the A/ = f selection rule. 

The same approximations in Eq. (ED, when applied to e, lead to 


e ~ e + i 


. IrriT c 
ReA 


(71) 

(72) 
-> 2 ir 

(73) 


Let us next go back to the mass matrix equations in (^1[) which, expanding in 
powers of e, we can rewrite as follows 


ReMi2 — ^ReLi2 

2 e ~ 


.Im M 12 - flrrir 12 


(74) 


i(Am+|Ar) |(Am + |Ar) ' 

To a first approximation, neglecting CP-violation effects altogether, we find that 

Am , „ „ Ar 


ReMi 2 — 


and 


ReT 12 ~ 


(75) 


If furthermore, we restrict the sum over intermediate states in r ]2 [see Eq. (pA|)1 
to 27r-states, an approximation which we have already seen to be rather good [see 
Eqs. (|T7p and (|35|)] we can write 

T 12 ~ (-iA* 0 e iS °)*iA 0 e iS ° = -(ReA 0 + iH) 2 , (76) 


from where it follows that 


ImT 12 2ReAoI m ^o „ hrn4 0 
Rcr 12 Re^o + IiilAq Rg-Aq 


(77) 


Then, using the empirical fact that Am ~ and U < Tg, we finally arrive at 
the simplified expression 


1 /. ImM 12 ImAo \ 
1 + i V Am Re A 0 / 


(78) 
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and, using Eq. (|73[), 


(79) 


1 - 7r /ImMi2 ImAo\ 

—e 4 -- -)- 

V Am ReA 0 / 

This is as much as one can do, within a strict phenomenological analysis of the 
CP-violation in K decays. We have reduced the problem to the knowledge of two 
parameters: e in Eq.(:7D|), and e' in Eq.(jT^). We shall come back to these parameters 
in Sec. 5. There, we shall discuss what predictions for these fundamental param¬ 
eters can be made at present within the framework of the Standard Model. As 
we shall see, the main difficulty comes from the lack of quantitative understanding 
of the low-energy sector of the strong interactions. In terms of QCD, the sector 
in question is the one of the interactions between the states with lowest masses: 
the octet of the pseudoscalar particles (n,K,rj). It seems therefore appropriate to 
examine the possibility of describing the interactions of these particles within the 
framework of an effective Lagrangian of the Standard Model at very low-energies. 
This will be the subject of the next three lectures. 


2 Introduction to Chiral Perturbation Theory 

Chiral perturbation theory (yPT) is the effective field theory of quantum chromo¬ 
dynamics (QCD) at low energies. In this lecture I shall give an overview of the 
basic ideas behind the yPT-approach to hadron dynamics at low energies. This 
will give us the basis for a study of possible applications of the yPT-approach to 
the non-leptonic weak interactions of K particles, which will then be the subject 
of the next lecture. 

There have already been several sets of TASI-lectures in previous years dedicated 
to ,\PT. The reader should consult them for complementary information lI 13. Some 
recent review articles can be found in | andjnj 


2.1 An Overview of the Basic Ideas in the %PT-Approach. 

In the limit where the heavy quark fields t, b, and c are integrated out, and the 
masses of the light quarks u, d and s are set to zero, the QCD Lagrangian is 
invariant under global SU(3)l x SU(3) r rotations ( Vl,V r ) of the left-handed and 
right-handed quark triplets 

q L = —-—q and q R = —-—q ; q = u, d, s. (80) 

> V L q L and q R —> Vr q R . (81) 

Formally, the global symmetry of the QCD-Lagrangian is in fact larger. The full 
symmetry group is SU(3) R x SU[3) R x U(l)v x U(1)a- The U(1)a symmetry is 
broken however at the quantum level by the abelian axial anomaly; the Adler 0- 
Bell-Jackiw 111 anomaly. The U(l)y quark-number symmetry is trivially realized 
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in the mesonic sector. At the level of the hadronic spectrum, this chiral 517(3) 
symmetry of the QCD-Lagrangian is however not apparent. Although the low-lying 
hadronic states can indeed be neatly classified in irreducible representations of the 
famous £77(3) symmetry of the Eightfold Way H-l, there do not appear degenerate 
multiplets with opposite parity. In QCD it is therefore expected, (and there are 
some good theoretical reasons for it H-j; as well as numerical evidence from lattice 
QCD simulations IL3) that the chiral-SW(3) global symmetry is spontaneously broken 
down to the diagonal SU( 3 )v=l+r group of the Eightfold Way. This pattern of 
spontaneously broken symmetry implies specific constraints on the dynamics of the 
strong interactions between the low-lying pseudoscalar states ( 7 r, K. ?/), which are 
the massless Nambu-Goldstone bosons associated to the “broken” chiral generators. 
As a result of the spontaneous symmetry breaking, there appears a mass-gap in 
the hadronic spectrum between the ground state of the octet of 0“-pseudoscalars 
and the lowest hadronic states which become massive in the chiral limit m u = rria = 
m s = 0 ; i.e., the octet of 1“ vector-meson states and the octet of 1 + axial-vector- 
meson states. 

The basic idea of the yPT-approach is that in order to describe the physics 
at energies within this gap region, it may be more useful to formulate the strong 
interactions of the low-lying pseudoscalar particles in terms of an effective low- 
energy Lagrangian of QCD, with the octet of Nambu-Goldstone fields ( A are the 
eight 3x3 Gell-Mann matrices, with trA a A fe = 2 S ab ) 


$0*0 = ■ ~${x) = 


( 7 r °/\/2 + 77 / a /6 


7T 

K~ 


7 r + K + \ 

-7r° / \/2 + 77 / \/6 K° 

K ° —277/ a /6 ) 


(82) 


as explicit degrees of freedom, rather than in terms of the quark and gluon fields of 
the usual QCD Lagrangian. In the conventional formulation, the Nambu-Goldstone 
fields are collected in a unitary 3x3 matrix U(x) with det U = 1, which under 
chiral-S'I/(3) transformations (Vl,Vr) is chosen to transform linearly: 


U - V R UVl 


(83) 


The effective Lagrangian we look for has to be then a sum of chirally invariant 
terms with increasing number of derivatives of U. For example, to lowest order in 
the number of derivatives, only one independent term can be constructed which is 
invariant under (IQ. Vr) transformations : 

4s = (84) 

where, as I shall soon explain in more detail, the normalization f n is fixed in such 
a way that the corresponding axial-current with the quantum numbers of the pion, 
obtained from this Lagrangian, is the one which induces the 7r —> pu transition. An 


12 




explicit representation of U is 


U (x) = exp 



A • 7ft(x) 


(85) 


and, from experiment, 

U - 92.4MeV. (86) 

Because of the non-linearity in <p, processes with different numbers of pseudoscalar 
mesons are then related. These are the successful current-algebra relations lid of the 
60’s which the effective Lagrangian formulation above EH2I incorporates in a compact 
way. 

The low-energy effective Lagrangian describes physical amplitudes by means 
of the lowest order term in a Taylor expansion in powers of momenta: 

A(pi,P2,.~) = J2 a ij ) Pi'Pj + °(P 4 )- ( 87 ) 


The expansion has no constant O(jP) term, due to the fact that the lowest order 
Lagrangian in (|84j) has two derivatives. As an exercise you should check some 


examples, like elastic 7r + 7r u scattering, where the induced amplitude one finds is: 


•4(P-hPO,P+,Po) = 


(p+ - p+y 

f 2 

J 7 r 


( 88 ) 


We expect this effective Lagrangian description to be useful for values of invari¬ 
ants of momenta sufficiently small as compared to the scale A x where spontaneous 
chiral symmetry breaking (SySB) occurs in QCD: 


p 2 /A 2 x < 1. (89) 

It seems also reasonable to expect A x to be of the same order of magnitude as the 
masses of the lowest states which become massive due to SySB, i.e., 


M p (770MeV) < A x < M A (1260MeV). (90) 

As we shall see, these expected features are justified phenomenologically. 

It is convenient to promote the global chiral SU (3) symmetry to a local SU (3) l x 
SU(3)r gauge symmetry. This can be accomplished by adding appropriate quark 
bilinear couplings with external field sources to £qc D , the usual QCD Lagrangian 
with massless quarks: 

Cqcd{x) = £q } CD ( x ) + ^(x)7 M lvp(x) +j 5 a^(x)]q(x) - g(x)[s(x) - ij 5 p(x)]q(x). (91) 

The external field sources v ft] a ;/ , s and p are Hermitian 3x3 matrices in flavour, 
and colour singlets. Under chiral-S'I/(3) gauge transformations (Vl(x), Vr(x)), they 
are required to transform as follows: 

Q = Vp -ap V L lpVl + iV L dpVl , 
rp = Vp + ap ->• V R rpV ] R + iV R dpV^ 
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(92) 

(93) 




and 

s + ip —> V R (s + ip)Vl- (94) 

( 2 ) 

In the presence of these external field sources, the possible terms in L vtl with the 
lowest chiral dimension, i.e., 0 (p 2 ) are 

4. = \fl {tr D„UDfU' + tr ( x U f + U X ')}, (95) 

where denotes the covariant derivative 

D,,y = df,U - ir^U + iUl (96) 

and 

X = 2 B[s(x)+ip(x)\, (97) 

with B a constant, which like f~. is not fixed by symmetry requirements alone. 
Once special directions in flavour space (like the ones selected by the electroweak 
Standard Model couplings) are fixed for the external fields, the chiral symmetry is 
then explicitly broken. In particular, the choice 

s + ip = M. = diag ( m u , m* m s ) (98) 


takes into account the explicit breaking due to the quark masses in the underlying 
QCD Lagrangian. 

Notice that the quark bilinear (f L {x)q l I j(x) has the same transformation prop¬ 
erties, under chiral-S'I/(3), as the U (x)-matrix field. In fact, U(x) can be viewed 
as the low-energy parametrization of the Nambu-Goldstone field excitations of the 
vacuum. To lowest order in the chiral expansion, the constant B appears then as 
the parameter which fixes the relative normalization to the light quark condensate: 

(0|gV|0) -fiBSij. (99) 


The relation between the physical pseudoscalar masses and the quark masses, 
to lowest order in the chiral expansion, is fixed by identifying quadratic terms in y 


in the expansion of the second term in (^5[), with the result 


X = 


m; 


+ M 2 


K+ 

0 

0 


M 2 K o 


mt 


0 

-M 2 

0 


K+ + M K° 


-mz 


0 

0 

+ Mf 


K+ 


+ M 2 


( 100 ) 


K° . 


The mass matrix in the 7 r°, 77 basis is not diagonal, because of a small admixture 
between these two fields. The appropriate diagonalization leads to the mass eigen¬ 
values: 


M 2 0 = 2 mB-e + 0{e 2 ), 

M 2 = ^-{fh + 2m s )B + e + 0(e 2 ), 
o 
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where 


and 


e = 


B (m u - m d f 


m x 


771 ) 


m = -(m„ + m d ). 


(103) 


(104) 


There are a number of important relations which follow from these lowest order 
results: 


The Gell-Mann-Oakes-Renner relation h—ll 

= —771 < 0 |uu + dd\0 > . 

The Current Algebra mass ratios B § 

M k + 


B = 


Ml 


Ml 


3 K 


2 irk (m u + m d ) (m d + m s 

The Gell-Mann-Okubo mass formula 0’ 0 


(2 TTi + 4m s ) 


3 Ml = AMI ~ Ml 


(105) 


(106) 


(107) 


Barring the possibility that the parameter B may be unexpectedly small (see 
however lj. and references thereof,) and neglecting the small O(e) effects, we are 
led to two quark mass ratios estimates 

Ml o - Ml + - (Ml -MI + ) em = 0 2g _ (iog) 


m d - m* 


m d + m h 


and 


m s — m 
2 rh 


Ml 

7 r u 

Ml 0 - Ml 


M 2 

7T0 


= 12 . 6 , 


(109) 


where we have subtracted the pion squared mass difference, to take into account 
the electromagnetic contribution to the charged pseudoscalar mesons self-energies. 
This electromagnetic contribution, in the chiral limit (' m U)d)S = 0,) gives a common 
mass to the charged K + and 7r + mesons. 

As already pointed out, in the case of the Standard Model, not all the external 
gauge held sources l and correspond to physical gauge fields. In the same way 
that the scalar and pseudoscalar sources are frozen to the quark mass matrix, as 
indicated in (|98|), the explicit chiral symmetry breaking induced by the electroweak 
currents of the Standard Model corresponds to the following choice: 


= eQ R [Afj, - tan i) w Z tl \; 


l = eQnlA^ - tan '0 W Z fl ] + 
e 


-Q 


( 3 ) 7 

L 


+ 


y/2 sin d w 


sin idw 

Bi +) wM + Q l (-)W<-)]. 


( 110 ) 


( 111 ) 
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Here the Q’s are the electroweak matrices: 


Ql = Qr = Q = -diag(2, -1, - 

1 ), 

Q { L = diag(l, —1, —1); 

(112) 

q£ = (eirb = 

( 0 

tud bus \ 


0 

0 0 . 

(113) 



0 0 / 



We can now compute, in a straightforward way, the chiral realization -to lowest 
order in the chiral expansion- of the electroweak currents of the Standard Model. 
They follow by taking appropriate variations of the lowest order effective action 
r( 2 ) = /d 4 x£W(x), with respect to the external field sources: 

A/h 2 ) j 

11 = VtlL = — = T 2 fl(D,U')U , (114) 

xr ( 2 ) j 

J“r = IrYir = -^- = 2 fl( D R u ) u ' • < 118 ) 

Expanding U in powers of $-matrix fields [see Eqs.(|82|) and (|S5|),] we have: 

Jt = + 0 ($ 3 ), (116) 

J £ = - -[$(1^$) - + C>(<h 3 ). (117) 

v 2 2 

Let us now consider the decay 7 r + —> // + z / /7 we mentioned earlier as an example. 
The hadronic matrix element can be readily computed using the form of the axial 
current deduced from the expressions just above, with the result 

< 0 |(J^)i 2 k + >=iV2f n p^ (118) 

explicitly showing the physical meaning of the f n constant in the chiral Lagrangian. 

You are now in the position of being able to calculate any observable you wish 
to lowest order in the chiral expansion; but I am sure you are more ambitious. How 
does one calculate higher order chiral corrections? This is what we shall now discuss 
in the next subsection. 

2.2 Chiral Perturbation Theory to 0(p 4 ) 

In QCD, the generating functional Z[v,a, s,p\ of the Green’s functions of colour 
singlet quark currents, is defined via the path-integral formula 

exp {iZ[v, a, s,p]} = J VqVqVG^ex p ji J d A xC QC D j , (119) 
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with jCqcd defined as in eq.(|9l[). The physical Green’s functions of a specific flavour 
are then obtained by functional derivatives with respect to the appropriate external 
field sources v, a, s, and p. The chiral symmetry properties which we have discussed 
imply that, at sufficiently small energies, there exists an effective Lagrangian £ e s of 
the Nambu-Goldstone field modes alone, in the presence of external field sources, 
such that 

exp {iZ[v,a, s,p]} — J VU[<&(x)\ exp ji J d 4 x£ c fr [G; v, a, s,p]|. (120) 

To lowest order in the chiral expansion, O(jr) as we have seen, the generating 
functional reduces to the classical action 


Z {2) [v,a,s,p] = J #xd%(U-,v,a, s,p), 


( 121 ) 


with as defined in eq.(|95l). 

To next-to-leading order in the chiral expansion, i ,e.O(p 4 ), the generating func¬ 
tional Z[v, a, s, p\ gets contributions from three different sources: 

• The most general local effective chiral Lagrangian of O(p 4 ). 

( 2 ) 

• The one-loop functional generated from the lowest-order £ ( , fl - Lagrangian. 


The well known Wess-Zumino-Witten functional induced by the 
abelian chiral anomaly El. 


non 


2.2.1 The Chiral Lagrangian of 0(p A ) 

The ingredients we have at our disposal to build this Lagrangian, their chiral trans¬ 
formation properties, as well as their chiral power counting are as follows: 


U, 

v R uvl , 

O(p 0 ) 

D„U, 

VrD.uvI , 

0 (p ); 

X, 

VrxVL 

o{p 2 ) 


VlF^vI , 

o(p 2 ) 


VrF^vI , 

o(p 2 ) 


Here, F'f"' and F'/T are the strength field tensors associated to the external 1^ and 
d u F-i[F,F], and Fj^ — <9 / V !/ — — i[r M ,r"]. (122) 


sources: 


F£ v = dn v 


Furthermore, since we shall only use at tree level, we can eliminate possible 
terms using the 0 (p 2 ) equations of motion obeyed by the [/-functional: 


(□£/)[/» - U(nu') = xU * - Ux 1 - fAxU' - U X ') 


(123) 
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Remember also that traceless 3x3 matrices Ai obey the constraint 

trAiA 2 A 3 A 4 A trA 4 A 3 A 2 A 4 A trA 4 A 4 A 3 A 2 
+trA 4 A 2 A 4 A 3 + trA 4 A 3 A 4 A 2 + trA 4 A 4 A 2 A 3 + 

—trA 4 A 2 trA 3 A 4 — trA 4 A 3 trA 2 A 4 — trA 4 A 4 trA 2 A 3 = 0. (124) 

Other constraints which we must also impose are that only terms which are in¬ 
variant under parity and charge conjugation should be allowed^ The most general 
Lagrangian which satisfies all these conditions is the following E3 

= U tr (D^D^U ) 2 + L 2 trD^DMtrD^D^U + 

L 3 trD^U^ D^UD u U^ D u U + 


L 4 tiD^D^Uti (x ] U + U^x) + L 5 tr D^WD^U (pU + U^ X ) + 


L 6 [tr (\ ] U + ^y)] 2 + L 7 [tr (pU - U ] x) 
L& tr (Ux'U X ' + rfxtfx) + 


+ 


iL 9 tr [F^D^UD V U ] + F^D^D U U) + L 10 tr U^Fj^UF Llu/ + 

Hi tr {F^F n , w + F^F J41V ) + H 2 tryty. (125) 

In this Lagrangian the parameters L^, i = 1,2,3,..., 10 are dimensionless coupling 
constants, which like f n and B in the lowest order effective Lagrangian, are not fixed 
by chiral symmetry requirements alone. The terms proportional to the coupling 
constants F{\ and H 2 involve only the external fields. As a result these coupling 
constants cannot be fixed from low-energy observables alone. Only if we had a 
detailed knowledge of the dynamics of how the effective chiral Lagrangian emerges 
from the underlying QCD Lagrangian could we give an unambigous definition of 
H i and H 2 . By contrast, as we shall later discuss, most of the other couplings 
can be fixed from low energy observables. Of course, in QCD, the L, constants, 
much the same as f n and B. are in principle calculable parameters in terms of the 
intrinsic A qcd scale only. We shall come back to this important question later in 
Sec. 4. For the time being we shall only be interested in the phenomenological 
determination of these constants. However, for that, we need to discuss first the 
possible contributions to low-energy observables from chiral loops. 


2.2.2 Chiral Loops 

Simple power counting tells us that loops generated by the lowest order Lagrangian 
are badly divergent. This is not a surprise: the non-linear sigma model in 4- 
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dimensions is not renormalizable i.e., an infinite number of local counter-terms are 
required. Here, however, we are considering the chiral Lagrangian as an effective 
field theory for low energies. Order by order in the momentum expansion of the 
effective theory, it is possible to specify a renormalizable framework. Of course, as 
we go to higher and higher powers of momentum, more and more local countert¬ 
erms will appear, with couplings which are not fixed by chiral symmetry properties 
alone, and eventually the predictive power of the effective field theory formulation 
will therefore disappear. To define the loop integrals it is necessary to fix a reg¬ 
ularization which preserves the symmetries of the Lagrangian. The well known 
dimensional regularization technique does the job beautifully. Since by construc¬ 
tion, the 0(p 4 ) Lagrangian contains all possible terms which are allowed by 
chiral invariance, all the one loop divergences -which by power counting can only 
give rise to local 0 (p A ) terms- can all be absorbed by suitable renormalizations of 
the Li and H l 2 constants. This program has been explicitly realized in the papers of 
Gasser and Leutwyler §1011. In fact, these authors have explicitly constructed the 
one-loop functional z| 0 4 ] [u, a, s,p\ at the required level of non-locality needed for 
all practical calculations. They do the path integral around the classical functional 
f/(<f> c i[u, a, s,p\), solution of the equations of motion in (|123|) with the boundary 
condition 1/(0) = 1. The method consists in expanding £ ef {(U-,v,a,s,p ) around 
$ = <f> c i and then doing the functional integral over the fluctuations £ = ( l> — $ c i. 
In order to obtain the one-loop effective action it is sufficient to expand to 

O[(0 2 ]: 

rS = / A4ff[4>ci] - f K, Di) + OK 3 ). (126) 

where D is a known, but rather complicated, differential operator acting on the 
fluctuating ^-variables. The one-loop functional is then embodied in the master 
formula: 

4 4 o P K a , s, p] = ^tr log D. (127) 

It is relatively easy to extract the singular part of Z^l p [v, a, s,p\. It corresponds 
to a local effective Lagrangian, exactly like the one in eq.( |125|) , with couplings L t 
and Hj\ 

L l i° p = Fi A loop , i — 1,2, 3, ...10; H] oop = f , A loop , j = 1, 2, (128) 

where 

A '“"‘ = £{dVI“^ |l0E(4,r)+r,(1) + 11 }’ 2 = 1.2; (129) 


and Tj, h y have the following rational values : 


r, = 


3 

0 ’ 


r 2 = 


3 

16’ 


f"3 — 0 , 


r 4 = 


i 

8 ’ 
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r 5 = 
r 9 = 


3 

8 ’ 

1 

4 ’ 



r 7 — o, 



r 8 = 


r 2 = 


5 

48’ 

5 

24’ 


( 130 ) 


The renormalized couplings L\(p) and //)(//), are then defined by subtracting from 
the tree level L t and Hj the one-loop singular contributions: 


U = L\{n) + Vi A Ioop , and M, = H](n) + f, A loop . (131) 


The renormalized coupling constants depend of course on the scale /i introduced by 
the dimensional regularization. The running in /i is governed by the coefficients T, 
(and Vj), which play the role of one-loop /^-functions: 


L'M = L'(/P + 4b logh (132) 

The fi -scale dependence cancels however in the full 0([V) calculation of a given 
physical observable. The non-polynomial contribution to a specific physical pro¬ 
cess will in general have a logarithmic p-scale dependence -the so called chiral 
logarithms- which cancels with the /j dependence of the tree level contribution 
modulated by the Lj(/x)-constants. Let us consider a typical example to illustrate 
this feature: the electromagnetic mean squared radius of the pion. The pion electro¬ 
magnetic form factor, when expanded in a Taylor series in powers of the momentum 
transfer, has the following structure: 

F 7T ±(q 2 ) = 1 + \ < r 2 >„± q 2 -\ -. (133) 

6 

The result of the calculation of < r 2 > n ±, to lowest non-trivial 0(p 4 ) in yPT is @ 


< r 2 > T ±= 12 


ft 


1 n M l 1 , 

16t r 2 / 2 ^ ° g p 2 + 2 ° g 


Ml 

A 2 


+ ? 1 . 


(134) 


The first term in the r.h.s. shows the tree level contribution from the renormal¬ 
ized Lagrangian; the rest of the terms are generated by the -finite though 
// dependent- contribution from the Feynman diagram loops (one-pion loop and 
one-kaon loop,) generated by the lowest Lagrangian. The scaling law of the 
L r 9 (^-constant, as shown in ( |132| ), is the same as the one of the chiral logarithm, 
and therefore the whole contribution is scale invariant, as it should be. 

There are a number of interesting generic features which emerge from this ex¬ 
ample, and which we next wish to point out: 


2 

• The factor is a characteristic factor of the loop-expansion. More pre¬ 

cisely, this factor appears modulated by the number of active flavour loops 7 if. 

2 

Therefore, we expect chiral loops to contribute O [rif x log’s] to physical 
processes in general. 
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• Experimentally 0 < r 2 > 7T ±= (0.439±0.008)fm 2 , showing that the contribution 
of the LI, constant, for any reasonable value of /i, dominates the electromag¬ 
netic mean squared radius of the pion. In other words, in this example, the 
tree level 0 (p 4 ) contribution largely dominates the “chiral logs” induced from 
the lowest order Lagrangian. This is a fact which can be easily understood 
within the framework of the 1 /iV c -expansioiT : in the large-JV C limit, Lg and 
/ 2 are 0(N C ); therefore the chiral loop contribution is 1 /Ay-suppressed with 
respect to the tree level contribution. 

• The phenomenological value we get, from this observable, for the Lg coupling 
constant at the scale of the p-meson mass: Lg(770MeV) = (6.9 ±0.7) x 10 -3 is 
in the range of the expected order of magnitude if, as we have assumed a priori, 
the chiral Lagrangian is a good effective theory for energies below the SySB 
scale A x . With a factor I/ 2 pulled out from , (the same -J'z factor which 
modulates the lowest order T) :fr Lagrangian,) the dimensionless couplings L ,, 

, I f2 

can then be traded by couplings, such that L\ = • It seems reasonable 

i x i 

to expect the A, : - constants to be of the same order of magnitude as A x . Then, 
from (|90|), we expect Af p (770MeV) < A, < My 1 (1260MeV); i.e., L\ ~ 10 -3 to 
5 x 10~ 3 . We see that L r g falls in this expected range; and in fact, as we shall 
soon see, so do all the other Lj-couplings. 

There is a wealth of experimental information on low energy hadron physics, 
which has permitted the phenomenological determination of the L, constants. The 
power of the chiral approach is that once these constants have been fixed from some 
experiment, or spectral function sum rules, there are of course predictions -and 
therefore tests- for other observables. In Table 1, I have collected the most recent 
compilation of the Li s. One can also read in the same Table 1, the experimental 
source which has been used for the determination of the appropriate constant. It 
would be nice to improve the accuracy of some of the low-energy experiments. 
Hopefully, the DATNE project at Frascati will eventually provide some of this 
improvement. A lot of theoretical work has recently been made in view of this 
project. This work has been published as a reference guide @. One can find there 
un update of many recent phenomenological applications of yPT. 

2.2.3 The Non-Abelian Chiral Anomaly 

Although the QCD Lagrangian with external sources is formally invariant under 
local chiral transformations, this is no longer true for the associated generating 
functional. The anomalies of the fermionic determinant break chiral symmetry at 
the quantum level. The anomalous change of the generating functional under an 

* This is the approximation proposed by ’t Hooft 0, where the number of colours N c in QCD is 
large with the product a s N c kept fixed. 
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Table 1: Phenomenological values of the renormalized couplings L\{M P ). 


i 

Li(Mp) x 10 3 

Source 

1 

0.7 ±0.5 

/P e 4 , 7T7T — > 7T7T 

2 

1.2 ±0.4 

A e4 , 7T7T — > 7T7T 

3 

-3.6 ±1.3 

A e4 , 7T7T — > 7T7T 

4 

-0.3 ±0.5 

Zweig rule 

5 

1.4 ±0.5 

Fk ■ F n 

6 

-0.2 ±0.3 

Zweig rule 

7 

-0.4 ±0.2 

Gell-Mann-Okubo, L 5 , L$, Sum Rules 

8 

0.9 ±0.3 

M k 0 - M k +, A 5 , (m s - rh) : {m d - m u ) 

9 

6.9 ±0.7 

lr 2 Y 

V /em 

10 

-5.5 ±0.7 

Ti — > evp 


infinitesimal chiral transformation 



V LjR = 1 ± ia T i/3 + ... 

(135) 

is given by E 

a 

5Z[v,a,s,p\ = ^ J 

d 4 x tr/3(x) fl(x), 

(136) 

where 

Q(x) = e pvap [v pu v ap ± ^ V p a 

2 

c jCLp ~i \vdfjCip } 



8 . 4 

“1”' 7^ fil'dp + — Qj^Qj U 

■OffOp], £0123 — 1; 

(137) 

and 

v pu = d p v v - d v v p - i [v p , v„], 

^^^ \VH ? Q'vl ’ 

(138) 


This anomalous variation of Z is an 0(p 4 ) effect in the chiral counting. 

Chiral Symmetry is the basic requirement to construct the effective yPT La- 
grangian. Since chiral symmetry is explicitly violated by the anomaly at the fun¬ 
damental QCD level, one is forced to add an effective functional with the property 
that its change under chiral gauge transformations reproduces ( [1361) . Such a func¬ 
tional was first constructed by Wess and Zumino cil. An interesting topological 
interpretation was later found by Witten 0. The functional in question, has the 
following explicit form: 

m t, r] wzw = - ±± j tr {EfE^E^EfE,^} 

- ^/ d‘xe„ m e [w(U,i,rY-“> - W(l,t,rY'°e) , (139) 
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with 


W{U,£,r)^ aP = tr {u£/J a tfr p +^U^r„Ue a tfr p + iUd/J a tfr p 

+ id^Utatfrp - i^e v U^r a U£p + E L ^d v r a Ul p 

- z L ^tfr a m p + zfadjp + xfaejp 

- I tftjjf, + - iE^E^j 

- (L R), (140) 


where 


s; = u'dM, 


£? = UdM\ 


(Mi) 

and (L <-> R ) stands for the interchanges U <-> W, <-»• and E^ <-> E^. The 
integration in the first term of Eq. ( |139| ) is over a five-dimensional manifold whose 
boundary is four-dimensional Minkowski space. The integrand is a surface term; 
therefore both the first and the second terms of T WZ w are O (p 4 ) according to the 
chiral counting rules. 

Since the effect of anomalies is completely calculable, their translation from 
the fundamental quark-gluon level to the effective chiral level is unaffected by 
hadronization problems. The anomalous action ( |139| ) has no free parameters. It 
is responsible for the 7r° —> 2y, 7 ] —> 2y decays, and the 737T, 77 t + 7t ~r] interactions 
among others. The five-dimensional surface term generates interactions among five 
or more Goldstone bosons. 

The variation of the anomalous action T[U, £, r) W zw with respect to appropriate 
external field sources, generates the chiral realization of the anomalous electroweak 
currents of the Standard Model, much the same as we calculated the electroweak 


currents in eqs.( |114|) , ( |115| ) induced from the lowest order chiral Lagrangian. With 


L„ = iU'D^U, 


(142) 


and Fj' v . the strength field tensors 


associated to l p and 77 , we have 


£r wzw 


N r 

c r llVpCT 

48vr 2 


iL„L„L„ + {/■;(;, + i rfF*U, L„ 



(143) 


This anomalous current is in fact defined, up to a chirally non-covariant polynomial 
in the external fields l. r. Only the covariant anomalous current above is measur¬ 
able. Of particular interest for semileptonic decays of K mesons, is the strangeness 
changing current which, in the presence of external electromagnetic interactions, 
couples to the charged W ± in the Standard Model: 


= 

w anom. 


N r 


£ pupa tr Q 


± 


D , 


\/2 sin 9w 487 T 2 
,U ] D p UD a U ] U + ieF up \u ] D a U, Q + \^Qu} 


+ h.c., (144) 
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where, here, 


I),, - UJj 


(145) 


is the covariant derivative with respect to electromagnetism only, (Q is the charge 
matrix in ( |112|) ,) and F ja , = () jL A v — d v A ll the electromagnetic held strength tensor. 
The first term in the r.li.s. contributes to /d; 4 decays such as 


K + (p) -> 7T + (p + )7T (p_)e + u e , q = p 


P+-P- 


(146) 


via the vector current, (H denotes the conventional phenomenological parametriza- 
tion of the vector form factor,) 


V“ = ~H(q 2 ) T7T e^q''(p + +p-Y(p+ -p-Y 


MJ 


with 


#(«) = 


N r . 


4 Ml 


Experimentally il 


16vr 2 /^ 3 


= -2.7. 




threshold 


= - 2 . 


±0. 


(147) 


(148) 


(149) 


The second term in ( 144 ) contributes to radiative semileptonic K -decays. An up¬ 
date of the rich phenomenology of these processes can be found in the DA4>NE 
nhvsics handbook b-1 


3 Weak Interactions and Chiral Perturbation The¬ 
ory 

In this lecture we shall study the weak interactions of Ji-particles within the frame¬ 
work of yPT. We shall be particularly concerned with processes induced by virtual 
W ± -exchange between quark currents, in the presence of the strong and electro¬ 
magnetic interaction. A central issue in the study of the non-leptonic weak decays 
of Ji-mesons, is the question of the origin of the A I = 1/2 selection rule. Let me 
start by explaining what the problem is. 

There is experimental evidence, both from K decays and hyperon decays, that 
the rates of non-leptonic strangeness changing transitions AS = 1 with isospin 
change A/ = 1/2 are particularly enhanced. This phenomenological fact is referred 
to as “the A/ = 1/2 rule”. The explanation of this selection rule has been a con¬ 
tinuous challenge to theorists for the last three decades! There are good qualitative 
indications that in the Standard Model, the observed A I = 1/2 rule has a dynam¬ 
ical origin; but no clear quantitative description of the enhancement has yet been 
exhibited. 

To illustrate the problem let us try to make a simple guess of the strength of 
the transitions expected for K —> titt decays in the Standard Model. For such 
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processes the W-mass can be considered as infinitely heavy, and the exchange of 
the 11^--field between two charged quark currents, ignoring for the moment gluonic 
interactions between the two quark currents, is described by a simple effective four- 
quark Hamiltonian: 


n_/AS=l 


= ^V ud V^ s 4(s L 'y' J 'u L )(u L 'y fl d L ) + h.c. 


with e.g., 


(s l Yul) = ^s^ a \x)y 


- — u {a) 


(z), 


(150) 

(151) 


and where a denotes colour indices. Here Gf is the Fermi coupling constant 
(1.166 x 10 -5 GeV~ 2 ) and Wd, V(* are matrix elements of the Cabibbo-Kobayashi- 
Maskawa flavour mixing matrix, (see the lectures of Roberto Peccei 0.) At this 
level of approximation, the effective Hamiltonian appears as the factorized product 
of two quark currents. In the previous section, and to lowest order in the chiral 
expansion, we have explicitly worked out the realization of these currents in terms 
of pseudoscalar fields, [see Eqs.( |114| ) and ( |1 15| ) ,] with the result: 




~^2 Ud,,K+ 


2V2 l 

i 


2V2 1 


( tt °d„K + ) + V3 (r i d li K + ) + V2 (n+d^K 0 ) 

(K + dy) + V3 (K+d^rj) + V2 tt + ) 


+ •••; (152) 


{uiij,d L ) 


-—=f n d u TT- 
V2 M 


+ 


2 L 


V2 

V 

+ 

(K°d,K~) 

V2 

O 

1 

+ 

( K-d^K 0 ) 


+ ••• 


(153) 


We have all the ingredients now to calculate the K 
I = 0, 2 which we introduced in the Sec.l., [see Eqs. 


7T7r isospin amplitudes 
and ([51]),] with the result: 


A ° = ( m * - ■ < im ) 
A 2 = ~V M V^2f, [Ml - Ml) . (155) 


As already mentioned in Sec.l as well, experimentally, from the ratio of the decay 
rates r(Ag —> 7r + 7r - ) and r(Ag —> 7r°7r°), and with neglect of radiative corrections, 
one finds 



~ 22 . 2 , 

exp 


(156) 


25 





















i.e., a factor of sixteen larger than our factorization estimate! More precisely, taking 
also into account the experimental rate for K + —> 7r + 7T°, we are led to the conclusion 
that our first guess underestimates the A/ = 1/2 amplitude by a factor of eight, and 
overestimates the A/ = 3/2 amplitude by a factor of two. 

One may be tempted to conclude that the factorization assumption we have used 
to make our estimate must indeed be a very naive picture. It turns out however, 
that this is precisely the result predicted by the leading behaviour of the 1 /A/" c — 
expansion in QCD; i.e,. the limit where the number of colours N c is taken to be 
large with a s N c fixed till. We shall see later why the large Ay estimate fails in this 
case to give the right order of magnitude. 

3.1 Short—Distance Reduction to an Effective Four—Quark 
Hamiltonian 

In QCD, the composite operator 

Q 2 = 4(sL7 M «L)(hz/Y/A), (157) 

is not multiplicatively renormalizable. As first pointed out by Gaillard and Lee 0 
and by Altarelli and Maiani 0, the effect of gluon exchanges generates a new AS = 1 
operator (a and [3 are quark colour indices; summation over repeated Greek letters 
is understood) 

Qi = 4 , (158) 

which mixes with the standard Q 2 operator in ( | 1 57| ) under renormalization. It was 
later noticed by Vainshtein, Zakharov and Shifman Lt3, that the effect of the so- 
called “penguin” diagrams where one light quark line is attached by gluon exchange 
to a weak s — d self-energy like diagram , brings further AS = 1 operators [qn = 
|(l + 75 )q(x)}: 

Qs = 4 {s L Yd L ) , 

q=u,d,s 

Q t = iflYdi £ 

q=u,d,s 

Q 5 = 4 (s L y i d L ) J2 (7/i7,/ifi) , 

q=u,d,s 

Qe = Yj dfffp/lR , 

q=u,d,s 

which also mix with Q 2 and among themselves under renormalization. The opera¬ 
tors 

Q-=Q 2 ~Qi (163) 

and Qi , i = 3,4, 5, 6, induce pure A/ = / transitions, but only four of these opera¬ 
tors are independent, since 

Q- + Qs - Qi = 0 . (164) 


(159) 

(160) 
(161) 
(162) 
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Under SU(3)l x SU(3)r transformations, and the operators Qi, i = 3,4, 5,6, 
transform like (8^, 1#) operators, while the combination 


Q (27) = 2 Q 2 + 3Qi - Q 3 (165) 

transforms like a (27^, 1 _r) operator which induces both A/ = | and A/ = | 
transitions via its components 

« (27> = 5 (Of/? + 5 «f/?) , ( 166 ) 

where 


/0-T\ 

Qi/2 = (s L Yd L ) + {sL'fuL) {u L 'y, x d L ) 

+2 (s L yd L ) ( din^d L ) - 3 ( s L Yd L ) (sl7m s l) , (167) 


and 


Q3/2 = { s L^d L ) (ml 7 m«l) + (sl^ul) {u L ^^d L ) 

- {s L ^d L ) (d L , ~Y t _ l d L s J . (168) 

The operator Q (2,) is multiplicatively renormalizable and does not mix with the 
others. 

The procedure to go from the Lagrangian of the Standard Model to an effective 
electroweak Hamiltonian, where only the degrees of freedom of the light quark fields 
u, d, s appear 0, consists in using the asymptotic freedom property 01 of QCD to 
successively integrate out the fields with heavy masses down to scales p 2 < m 2 , i.e. 
below the charm quark mass. The appropriate technique is the operator product 
expansion bJ and the use of renormalization group equations Ldl to compute the 
various Wilson coefficient functions of the four-quark operators Qi, i = 1,..., 6. 
The inclusion of virtual electromagnetic interactions in the process of integrating 
out the fields with heavy masses brings in li-f H-J four new four-quark operators to 
the effective electroweak Hamiltonian. With e q denoting the corresponding quark 


charges in units of the electric charge, the new operators are: 

Qi = 6 {sL'fdL) Y e <i («R7/i9fl) > 

a=u.d.s 

(169) 

Qs 

= 6s£7 ^d ,3 L Y e q <&ln<lR, 

q=u,d,s 

(170) 

Q 9 

= 6 {SL'fdL) Y e Q feww) , 

q=u,d,s 

(171) 

Q 10 

= 6 s a L ^d! 3 L Y e v $TUAl ■ 

q=u,d,s 

(172) 
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Under the action of the chiral group SU(3)l x SU(3)r the operators Q 7 and Q$ 
transform like combinations of (8^, 1#) and (8^, 8r) operators, while Q$ and Q w 
transform like combinations of (8^, 1#) and (27 L , 1 R ). 

Owing to the big value of the top quark mass, higher-order electroweak contri¬ 
butions turn out to be relevant when analyzing some CP-violation effects. These 
additional corrections do not generate new operators, but contribute to the Wilson 
coefficient functions of the operators Qi, i — 1,..., 10. 

In processes with leptons in the final state, three more operators need still to be 
considered ffiJ: 


Q 11 

= 4 ( s L Yd L ) Y, , 

l=e,fi 

(173) 

Q 12 

= 4 {sLl^di) ^2 (IfCI^r) , 

(174) 

Q13 

= A{s L Yd L ) {vlIhVl) ■ 

(175) 


Since these operators are bilinear in the quark-fields, they have to be considered 
separately whenever they may be relevant. The operators Qn and Q 12 play an 
important role in the transition Kl —> 7r°e + e _ which we shall study in Sec.5. 

In the Standard Model, the AS = 1 non-leptonic interactions can then be 
described by an effective Hamiltonian 

s~i 10 

». f/ =1 = -mVM E C M Qi + h.c., (176) 

v 2 i=! 


where Qi are the ten local four-quark operators introduced above, and C,(ji) the 
modulating Wilson coefficients which are functions of the masses of the fields which 
have been integrated out, i.e. t, Z, W, b, and c, as well as of the overall renormal¬ 
ization scale /i. The Wilson coefficients of the ten four-quark operators take into 
account the effect of the strong interactions down to scales of O(n), in the presence 
of virtual electroweak interactions, which in practice are kept to O(a). 

The Wilson coefficients C',;(/i) obey the renormalization group equation — ck s ) 


d d 


C 




w 




,Oi s ,a) = 7 T (a s ,a)C 


Mi 


w 


2 1 




where P(a s ) denotes the QCD beta function 

PM = /3i— + /? 2 T—) 2 

7r \ 7r / 


+ 


Pi 

P2 


1(^11^+271/), 


— 17 fN c \ 2 1 Nl — 1 n-/ 

3 V 2 / 2 2 N c 2 
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H- N c n f . 

12 * 



(177) 


(178) 

(179) 

(180) 









The effect of electromagnetic corrections in the (3 function can be neglected at the 
level of precision that the Wilson coefficients are needed at present. The anomalous- 
dimension matrix 7 ^ is defined by the mixing of the operators (), under renormal¬ 
ization: 

( Qi) = ~\lij{a a , «) (Qj) • (181) 

In perturbation theory, the matrix elements of 7 which are needed for a calculation 
of the Wilson coefficients at the leading and next-to-leading logarithmic approxi¬ 
mation, are the coefficients of the first four terms in the expansion 


7 («s,«) = 


a ° Nc l „M _i_ « 1 (°) + 
7T 2 ls 7T2 7e 


a, N a 


7 r 


1 (i) , cn s N c a 1 

■jTs ^ TT TT A^ Se ^ ' ' ' ’ 

4 7T 7T 4 


(182) 


The most recent calculations of the full 10 x 10 anomalous dimension matrix corre¬ 
sponding to this expansion, can be found in Refs. 

The general solution of the renormalization group equation above is given by 

/ A/w - \ [ roisiiJ.) r y'i'( z ) -■ 

C = lT as ex p/ dz 1 ’ > C(l, a s (M w ), a), (183) 

V d J { Ja s {M w ) zf3{z) ) 


where T aa denotes ordering in the QCD coupling constant increasing from right to 
left. The vector C(l, a s (M w ), a), which defines the initial boundary conditions, has 
to be extracted from the calculation of the full Feynman diagrams at the M\y mass 
scale. In the limit a s = a = 0: C 2 (Mw) = 1 and all the other C\ vanish. To 0(a), 
the Ci with i = 1, 2, 3, 7 and 9 get contributions, while C 8 = Ci 0 = 0. The two-loop 
final expressions for the Wilson coefficient functions are rather elaborate. I shall 
not reproduce them here. They can be found in the references quoted earlier li-jLiJ. 

We can now understand why our factorization estimate,earlier, failed. The effect 
of gluon exchanges between quark current bilinears brings in, by mixing, new four- 
quark operators. The Wilson coefficients of the full set of four-quark operators has 
to be scaled from the large t mass and IF mass scales down to a renormalization 
mass scale /a ~ lGeV, below which the QCD perturbative evaluation is no longer 
trustworthy. This long running of scales brings in large logarithms, the effect of 
their size being governed by the anomalous dimension matrix of the four-quark 
operators. It is only when terms of at least O(N 0 ) are taken into account in the 
weak amplitudes, that the effect of the long evolution from short-distances to long- 
distances is incorporated. To see this more explicitly, it suffices to compute the 
lowest order Wilson coefficients in the sector of the Qi and Q 2 four-quark operators. 
In the basis of the Q- = Qi — Qo and Q + = Qi + Q 2 operators, you will find 


cV) 


/ ®s(Mw) \ 

V W(h 2 ) / 


±3/2(l T l/iVc) 

■4 


(184) 


The net effect is that transitions induced by the operator are enhanced by a 

/ m 2 2 \ F? 

factor log / log ~ri~ ) ; while the transitions induced by the Q + -operator are 
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depressed by a factor (log 


/ l MS, 
- / log — 21 


A ms 


y/7 


Although this does not explain the 


observed A/ = 1/2 enhancement, it goes well in the right direction, and points 
towards a solution of the problem: the possibility that the enhancement due to the 
anomalous dimensions continues in the evolution of the hadronic matrix elements at 
low energies. 

Another interesting issue, which appears at the two loop level of the perturbative 
calculations, is the question of renormalization scheme dependence. The Wilson 
coefficient functions have been calculated in two schemes: the ’t Hooft-Veltman 
renormalization scheme 0, and dimensional regularization with an anticommuting 
7 5 . The scheme dependence can only be removed by matching the calculation of 
the C/s with a similar calculation of the matrix elements of the Qi(n) operators 
in the same scheme. Unfortunately, the technology to calculate low energy matrix 
elements of light four-quark operators is not yet developed to the degree of sophis¬ 
tication of perturbative QCD. Only approximate methods have been developed so 


jrd Bl 


e the lectures of Steven Sharpe B) ; various versions 
the 1/N C approach developed by Bardeen, Buras 
and more recently, the QCD low-energy effective action ap- 
We postpone the discussion of some of these methods to the following 
sections, where we shall review estimates of various \PT coupling constants relevant 
to K decays. 


far: lattice models 
of QCD sumr 
and Ge. 
proach I 


3.2 The Effective Four—Quark Hamiltonian in %PT 

The effective four-quark Hamiltonian that we have derived in the previous sub¬ 
section contains only light-quark fields degrees of freedom. The light quark fields 
still have strong interactions mediated by the gluon fields of the QCD-Lagrangian. 
Now we want to find the effective chiral realization of this Hamiltonian in terms 
of (pseudo) Nambu-Goldstone degrees of freedom only; i.e., we want a description 
in terms of an effective chiral Lagrangian -with the same chiral transformation 
properties as the four-quark Hamiltonian- which will be incorporated as a weak 
perturbation to the strong chiral effective Lagrangian we discussed in Sec. 2. We 
also want to construct the weak effective Lagrangian in terms of a chiral expansion 
in powers of derivatives and quark masses, much the same as we did for the sector 
of the strong interactions. 

As we saw in the last subsection, there are essentially^ two types of four-quark 
operators, according to their transformation properties under chiral-S77(3): those 
which transform as (8/, > and those which transform as (27/,, 1#). To lowest 

order in powers of derivatives, all the possible operators one can construct with 
these transformation properties are the following: With C^x) the 3x3 flavour 

t There is also the combination_Qf the operators Q 7 and Q § which transforms like (8l,8r) which 
has to be considered separately c3. 
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matrix field 


(185) 


= -i^jU\x)D IL U(x), 
which by itself transforms as C f , —> VlC^VI, we can construct the two operators 

£ 8 (x) = E(^)»(^) a ; 


£-27(x) — g (£/J 2 l (^m)i3 + (AJ23 (AJll ' 


(186) 

(187) 


The second operator induces both A I = 1/2 and A I = 3/2 transitions via its 


components 


where 


0 


tn(x) = g£ 2 7 


<1/2) ' I )+“£! 3/21 ' 


-27 


X 


( 1 / 2 ), 
27 1 

(3/2) 
27 1 


X 


X 


~ (£( 1)21 (^ m ) i 3 + (^ m ) i 3 ^ 5 (£//) 22 

= (A ^2 :/ (^-m)i 3 (AJ23 (^)ll _ (AO22 


(188) 

(189) 

(190) 


The most general effective Lagrangian we are looking for, to lowest order in the 
chiral expansion, is then* 


r AS=l 

^eff 


» = 


gs£8 + S27^-t/ ^ + ~S27^27 
5 


( 1 / 2 ) 


+ h.c.. 


(191) 


The two constants g 8 and g 27 are dimensionless constants, real constants in as far as 
CP-violation effects are neglected, which like f n , B, and the L, coupling constants 
of the strong Lagrangian, cannot be determined from symmetry arguments alone. 
With the factors of f n included in the definition of C M in ( |185[ ), the constants g ^^7 
are of 0(1) in the large Ay expansion. We can get their phenomenological values 
from a comparison between the expressions for the K —> tttt isospin amplitudes Aj 
calculated with the Lagrangian above: 


/=' / u dKi(g 8 + L „)y/2f, (Ml - Ml) , 


A 0 — 

V 2 

Ai = ~K d V;,g 27 2 A (Ml - Ml) , 

and experiment; with the result 0: 


|gs 


g27 |exp. - 5.1, 


|g27 


exp. 


~ 0.16. 


(192) 

(193) 

(194) 


( A coupling ~ (xU { + UX r )'23 is also possible a priori. However, for on-shell processes, this term 
can be rotated away so as to maintain the normalization condition < U >= 1 . The physical effect is 
then of higher order. This term has still physical relevance, even to lowest order, for off-shell non- 
leptonic Green’s functions, and therefore brings in an extra coupling constant which is not fixed by 
symmetry requirements alone. 
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From a comparison between the ^/-amplitudes calculated above, and the results of 
our factorization estimate in Eqs.( |154|) and ( |155|) , we can also read the values of 
the coupling constants g 8j2 7 corresponding to this approximation: 

3 1 

§81 fact. = V, §271 fact. = ~■ (195) 

5 3 

rather far away, as already discussed, from the experimental values. 

Once the couplings gg and g 27 have been fixed phenomenologically, there follow 
a wealth of non-trivial predictions for K —> 777777 decays and some radiative K— 
decays as well. Concerning the latter, there appear some interesting features, which 
I think are worth pointing out, because they reveal the power and the simplicity of 
the chiral approach: 

• K-decay amplitudes with any number of real or virtual photons mid atrrwst one 
pion in the final state vanish to lowest 0(p 2 ) in the chiral expansion 

The reason for it is due to the fact that electromagnetic gauge invariance 
requires physical amplitudes to have a number of chiral powers higher than 
just the two powers allowed by the lowest order effective Lagrangian. Processes 
like 

K s 77 , K l -> 77°77; (196) 

K + —■> 7r + 77, Ks —•► 7r°77; (197) 

K + —■> 7r + e + e _ , Kl —> 7r°e + e _ , (198) 

are all of this type. They are at least O(j > 4 ) in the chiral expansion. 

• K-decay amplitudes with two pions and any number of real or virtual photons 
in the final state factorize, at 0(p 2 ), into the corresponding K —> 7777 on-shell 
amplitude times a universal bremfitrahlung-like amplitude. 

The proof consists in a simple adaptation of a well known theorem due to F. 
Low clll, to lowest order in ,\ PT. [A sketch of the proof for one photon is given 
in Ref. E3. Applications to K —> 77777 can be found in Ref. Lil] 

3.2.1 Weak Amplitudes to 0(p 4 ) in xPT 

The full analysis of the one loop divergences generated by the lowest order weak 
Lagrangian, in the presence of the strong effective chiral Lagrangian; as well as 
the classification of the possible local terms of 0(p 4 ) was first made in and, 
using different techniques in as well. The number of terms is too large to do a 
phenomenological determination of the couplings, as it has been done in the strong 
interaction sector. Even restricting the attention to the effective realization of the 
four-quark operators which transform like (8 l, It?), still leaves twenty two possible 
terms that contribute to non-leptonic K decays with possible external photons 
or virtual Z bosons. To determine phenomenologically these couplings is not a 
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question of calculation complexity; it is simply that there is not enough available 
experimental information to do the job! The only possible way to do yPT usefully in 
the sector of the non-leptonic weak interactions is to combine the chiral expansion 
with other approximation methods, like e.g., the l/lV c -expansion; or to resort to 
models of the low-energy effective action of QCD that one can first test in the 
strong interaction sector. 

The art of the game in making clean yPT predictions for non-leptonic weak 
processes, consists in finding subsets of observables which to 0(p 4 ) are fully given 


and Kt 


7T°77 


or which involve a 


For 


K —^ 7 rZ + Z decays 
I shall discuss some of these 


by a chiral loop only, like e.g., Kg —> 77 
small number of unknown 0(p 4 ) local couplings, li 
a recent review on the the state of the art see Ref. 
processes in Sec. 5. 

Another sector where it has been possible to test the validity of the chiral ex¬ 
pansion in non-leptonic weak interactions is in K —> 2i r and K —> 37 t decays. The 
description of these decays to 0(p 4 ) in yPT involves seven linear combinations of 
coupling constants. Altogether, imposing isospin and Bose symmetries, these decays 
can be parametrized in terms of twelve observables. Five of these parameters, the 
quadratic slopes in the Dalitz plot for the various K —> 3n modes, vanish to lowest 
order in the chiral expansion. The resulting five constraints can be formulated in 
terms of neat predictions lAl for the slopes. The five predictions are compatible with 

r important result which also emerges from the 
, and which is relevant for the understanding of 


experiment within errors. Anot 
K —> 27T, 37 t amplitude analysis 


the underlying dynamics of the A/ = 1/2 rule, is the fact that, in the presence 
of the 0(p 4 ) corrections, the fitted value for |g 8 | is ~ 30% smaller than the lowest 


order determination in Eq. 


This is, again, another little factor which, like the 


short-distance enhancement that we discussed earlier in subsec. 3.1, helps towards 
explaining the A/ = 1/2 rule, but the bulk of the enhancement remains still unrav¬ 
elled. Let me note that, at the same order of approximation, the coupling g 27 has 
110 large corrections. 


3.2.2 Penguins in the Large-N c Limit 


It has often been speculated that the bulk of the origin of the A I = 1/2 enhancement 
comes from the large matrix elements of the Penguin-generated Qq operator in ( |162|) . 
By Fierz reordering, this operator can also be written as 


Qq — = 8 55 ( s l<1r ) (< 7 /?^) • (199) 

q=u,d,s 

I find it very instructive to discuss the chiral realization of this operator El, because 
it reveals a number of interesting features; and also because it provides an excellent 
example of the way that, when systematically combined, the chiral expansion and 
the large Ay expansion may eventually help us to make substantial progress in the 
understanding of low-energy QCD. 
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Using Eqs.(|9lD and ©, the term which couples quark bilinears to external 
scalar and pseudoscalar sources can be written as follows 


- g(s - i-f 5 p)q = -7yg(qnX(iL + qLX ] qn)- 


( 200 ) 


In the large-JV C limit, the operator Qq in ( |199|) factorizes into the product of two qq- 
densities. The chiral realization of these densities to 0(p 2 ) in the chiral expansion 
and keeping only those terms which may eventually contribute to the lowest order 
gyU 8 - piece in the Lagrangian in ( |191| ), proceeds as follows: 


{s L qm) = 2 B—f 

°X3i 

(< qmd L ) = 2 B^ 
°Xi2 


23 + L 5 (UD^D^U) i3 + • • - j , (201) 

2D [xp U i + L b {D^D»UU') 2i + • • j . (202) 


Since UW = U^U = 1, there is no contribution to Q 6 of O(p 0 ), and to leading order, 
both in the chiral expansion and in the large-JV C expansion, we hnd 

Qe => -8 x 4 B 2 -L 2 L 5 (D^D,U) 23 . (203) 

^ J 7r 


We can cast this result in terms of the contribution to the gs-coupling, induced 
by the term in the effective four-quark weak Hamiltonian proportional to the Q§- 
operator: 


[S 8 ]Qta — Penguin ) 


< 'ifj'lp > 

J 7r 


(204) 


As I said before, there are a number of interesting features emerging from this result: 


• Remember that in the large- N c limit f 2 ,L§ and < vjvj > are all 0(N C ). 
Because of the fact that C§ ~ a s , [gs]6 th -Pengmn is next-to-leading in the 
large-IV c limit. It is precisely this next-to-leading contribution that we have 
succeeded in calculating explicitly. 

• The Wilson coefficient C& has an imaginary part induced by the CP-violation 
phase in the Cabibbo-Kobayashi-Maskawa matrix E3. It is precisely this con¬ 
tribution that, in the Standard Model, induces the e'-amplitude we discussed 
in Sec. 1. The fact that we can say something quantitative about the size of 
the modulating coupling [gs] 6 th -Penguin is of course welcome for phenomenology. 
We shall come back to this in Sec. 5. 

• In QCD the matrix element < > is /j -scale dependent. The scale depen¬ 

dence is given by the anomalous dimension of the '0'0-operator. On the other 
hand, the coupling constant g 8 is a scale independent quantity. This example 
exhibits explicitly the cancellation H between the p-scale dependence of the 
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short-distance Wilson coefficient Ce(^ 2 ) and the // scale dependence of the 
long-distance four-quark operator bosonization, which appears via the con¬ 
stant (< '0-0 >) 2 . In the large-IV c limit the operator Q e does not mix with the 
others because, in that limit, the anomalous dimension matrix in Eq.( |182| ): 

In that limit, C 6 depends on /i via C 6 ~ a(/i 2 ) 9/n , which 


7; 


(o) 


\NA 


*7 


C v 66 - 


exactly cancels the /i-dependance of (< vrS >) 2 evaluated in the same large- 


N c limit. 


• Of the parameters which appear in the r.h.s. of Eq.( |204| ), the one which has 
the largest uncertainty is < VJ'ip >. The most recent determination from an 
update of QCD sum rules methods gives the value E3: 

< -0-0 > (lGeV 2 ) = -(0.013 ± 0.003)GeV 3 . (205) 


It is also questionable which value of f n one should use in the r.h.s. of Eq.( p04|) . 
At the level of approximation that the calculation has been made, it seems 
appropriate to take the value of f n in the chiral limit; i.e., f.^ ~ 86 MeV. 
Then, using this value for f n , the valuers ~ 1.4 x 10 -3 which is in Table 1, 
and the result of the 6 / calculation in Ha, I find 

[Ss] 6 th — Penguin — 0.9. (206) 

From this result, one is led to conclude that, unless the next-to-leading or¬ 
der 1/IV C corrections to this calculation are huge, the bulk of the A I = 1/2 
enhancement does not come from the Penguin Q e operator. Barring this pos¬ 
sibility, the long-distance contribution from the Q_ operator remains then 
the most likely candidate to provide the bulk of the enhancement. 


3.3 AS = 2 Non-Leptonic Weak Interactions 

The Standard Model predicts strangeness changing transitions with AS = 2 via 
two virtual IP exchanges between quark lines, the so-called box diagrams. The 
reduction via the operator product expansion results in an effective Hamiltonian Eil 
which is proportional to the local four-quark operator 


Qas=2 = (sl7^l)(s7/^l), 

modulated by products of the flavour mixing matrix elements 


(207) 


\ = V* A V qs , q — u, c, t; (208) 

times functions Tig ,3 of the heavy masses 7n 2 ,M^,ml,m 2 of the fields which have 
been integrated out: 


n 


AS=2 

eff 


G 2 F Myy 
47r 2 


A 2 Fi + + 2A c At_F3 


Cas=2 


2 )Qas=2 ■ 


(209) 
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The operator Qas=2 is multiplicatively renormalizable and has an anomalous di¬ 
mension 7 (cn g ) defined by the equation 




dji- 


< QaS=2 >= -^7( a s) < QAS=2 > 


At the one-loop level 

jM = -7i + 0 (~) > 7i = |(1 ~ 1 / N c)] 

7T \ 7T / 2 


( 210 ) 


( 211 ) 


and 


c as , 2 (P) = 0i (PH - 31 ~ 2/9 


The matrix element 


< k°\q a ,= 2 (0)\k° >= - s\m\b k ( P), 


( 212 ) 

(213) 


defines the so-called Bk parameter, which governs K° — K° mixing at short dis¬ 
tances, and is one of the crucial unknown parameters in the phenomenological 
studies of CP-violation in the Standard Model. The definition above is such that 
in the so-called vacuum saturation approximation @ 


(B k )v s = 1. (214) 

In the large- Ay limit, the four quark operator Qas=2 factorizes into a product 
of the two left-handed current operators: (s^y^d^) times (s^y^d^). Each of these 
currents, to lowest order in yPT, has a bosonic realization as indicated in Eq. (|114| ); 
i.e., 

(s L y M d L ) =► H-. (215) 

r, x2 

This approximation results then in a value: B K = 773 -. There are however chiral 
corrections of 0(p 4 ) to this determination which are still leading in the large-lV c 
limit. The leading 1 /Ay-corrections come from the L r> term in as well as from 
the K 0 wave-function renormalization, with the result 

(sa“d L ) =c ^1 4 . (l - d„R -° + • • ■. (216) 

where the first factor in the r.h.s. is the one coming from . (Notice that 
the contribution from chiral logarithms is 1/N C suppressed.) The net result is a 
renormalization of f n . This renormalized f n coupling is a first approximation to 
the physical /#-coupling constant which is measured in K + —> /i + 7y, decays: fx — 
113MeV. To the extent that the two are identified, one then obtains the result 
which is often quoted in the literature as the large-N c prediction: 
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4 Models of the QCD Low-Energy Effective Ac¬ 
tion 

The purpose of this lecture is to give an overview on models of the low-energy 
hadronic interactions, which have been developed during the last few years, and 
which try to focus on general properties that the QCD low-energy effective action 
is expected to have. The aim here is to find generic features of these models which, 
eventually, one may be able to promote to the rank of a low-energy effective action 
derived from QCD, hopefully within a well defined set of approximations. 

Let me first formulate the problem which one would like to solve. For this 
purpose, it is convenient to promote the global chiral-Ff/(3) symmetry to a local 
gauge symmetry, in the same way that has already been discussed in Sec. 2. It is 
also convenient to use a path integral representation for the generating functional 
T (v,a,s,p) of the Green’s functions of quark currents: 

e iT(v,a,s,p) = 1 j VGij det fl exp (_ z J d 4 x 1 g^G^) , (218) 

with TJ) the Dirac operator 

P = + ig a Gp) - + 150 - 11 ) + - H5P ); (219) 

where G^ is the gluon field, G^ the gluon field strength tensor, and a M , s,p 
external field sources. The normalization factor Z is such that r(0, 0) = 1. 

The chiral symmetry of the underlying QCD theory implies that r(u,a, s,p) 
admits a low-energy representation 

e ir(« I o, SJ> ) = I j VU exp[i Jd*x £es(U ; v, a, s,p )], (220) 

in terms of an effective Lagrangian C e ^(U\ v, a, s, p) with U(x) a 3 x 3 unitary matrix, 
with det U = 1, which collects the octet of pseudoscalar fields ( 7 r, K, rj). 

There is only one term in T efr which is known from first principles. It is the 
term associated with the existence of anomalies in the fermionic determinant 
The corresponding effective action is the Wess and Zumino Lil functional that we 
have discussed in Sect. 2. All possible other terms in £ e ff, are not fixed by symmetry 
requirements alone. It would be nice to find, at least approximatively, a dynami¬ 
cal scheme rooted in QCD, with as few free parameters as possible (ideally A qcd 
only!), which allows for a derivation of the coupling constants in the effective chiral 
Lagrangian. As we have seen in Sec. 3, the need for such an approximate dynami¬ 
cal scheme is urgently needed at present to make progress in the phenomenology of 
non-leptonic flavour dynamics. 

The various types of models which have been discussed in the literature can be 
classified, roughly speaking, in one of the following entries: 

• QCD in the Large-iV c Limit. 
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• Low-Lying Resonances Dominance Models. 

• The Constituent Chiral Quark Model. 

• Effective Action Approach Models. 

• The Extended Narnbu and Jona-Lasinio Model ( ENJL-model.) 

4.1 QCD in the Large— N c limit 

It would be a major breaktrough, if one could derive the low-energy effective La- 
grangian of the interactions between Nambu-Goldstone modes in the large N c limit 
of QCD. So far, it has only been possible to obtain constraints among various cou¬ 
pling constants in this limit; but not their values in terms, say, of Aqcd- A typical 
example is the relation: 2 L\ = L 2 . which, as first noticed by Gasser and Leutwyler hll. 
follow in the large N c limit of QCD. Unfortunately, nobody can claim as yet to be 
able to compute , say L 2 , in that limit. Often in the literature, there appear state¬ 
ments about “ large Ay predictions” but, in fact, they have been all derived with 
some extra ad hoc assumptions. 

An interesting approach to do approximate calculations within the framework 
of the 1/Ay expansion is the one proposed by Bardeen, Buras and Gerard 000. 
which they have applied extensively to the calculation of non-leptonic weak matrix 
elements. The basic idea is to start with the factorized form of the four-quark 
operators in the effective weak Hamiltonian, and to do one-loop chiral perturbation 
theory, keeping track of the quadratic divergences which appear. If one was able to 
work with the full hadronic low-energy effective Lagrangian, it would be possible to 
obtain a smooth matching between the scale dependence of the Wilson coefficients, 
calculated at short-distances, and the hadronic matrix elements calculated with 
the full hadronic low-energy effective Lagrangian. The hope with the approach 
proposed by Bardeen, Buras and Gerard is that the numerical matching of the 
quadratic long-distance scale with the logarithmic short-distance scale, may turn 
out to be already a good first approximation to the problem one would like to solve. 
The technology of their approach is explained with detail in their papers. 

4.2 Low-Lying Resonances Dominance Models 

There has been quite a lot of progress during the last few years in understanding the 
role of resonances in yPT. At the phenomenological level 00. it turns out that the 
observed values of the Lj-constants are practically saturated by the contribution 
from the lowest resonance exchanges between the pseudoscalar particles; and par¬ 
ticularly by vector-exchange, whenever vector mesons can contribute. The specific 
form of an effective chiral invariant Lagrangian describing the couplings of vector 
and axial-vector particles to the (pseudo) Nambu-Goldstone modes is not uniquely 
fixed by chiral symmetry requirements alone. When the vector fields describing 
heavy vector particles are integrated out, different field theory descriptions may 
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lead to different predictions for the L t couplings. It has been shown however that 
if a few QCD short-distance constraints are imposed, the ambiguities of different 
formulations are then removed Ed. The most compact effective Lagrangian formu¬ 
lation, compatible with the short-distance constraints, has two free parameters: 
and M v . When the vector and axial-vector fields are integrated out, it leads to 
specific predictions for five of the L* constants: 

4 y) = 4 y) /2 = -4 y) /e = 44/8 = -44^4 = ^ o.e x 10 - 3 , ( 221 ) 

LbMy 

in good agreement, within errors, with experiment. [See Table 1.] 

It is fair to conclude that the old phenomenological concept of vector meson 
dominance (VMD) E3 can now be formulated in a way compatible with the chiral 
symmetry properties and the short-distance behaviour of QCD. 

In view of this success, there have been various suggestions of extensions of 
VMD-models in the literature. In particular, one would like to extend the idea of 
VMD to the non-leptonic sector of the weak interactions, in order to have a useful 
low-energy chiral effective Lagrangian formulation. (Remember the problem of the 
proliferation of couplings we have discussed in Sec. 3.) At the strict phenomenolog¬ 
ical level, one is forced to introduce weak couplings between vector (axial-vector) 
fields and/or the (pseudo) Nambu-Goldstone fields. The fact that these weak cou¬ 
pling constants are experimentally unknown, forces one to resort to models if the 
idea of VMD is to be pursued. 

Models like the Geometric Model E9E3 and the Quark Resonance Model EU are 
attempts to reduce the number of free coupling parameters which are allowed in 
principle by chiral symmetry requirements alone. However, the precise relation of 
these models to the specific assumptions which one is making within the underlying 
QCD-theory to justify them, remains unclear. 

4.3 The Constituent Chiral Quark Model 

This model was introduced by Georgi and Manohar E 2 . in an attempt to reconcile 
the successful features of the Constituent Quark Model 4 with the chiral symmetry 
requirements of QCD. The basic assumption of the model is the idea that between 
the scale of chiral symmetry breaking A x and the confinement scale rsj Aqcd the 
underlying QCD-theory, may admit a useful effective Lagrangian realization in 
terms of constituent quark fields Q ; pseudoscalar particles ; and, perhaps, “gluons”. 
The Lagrangian in question has the form 

^eff 1 = ’iQlnid/J. + + r 'fj)Q + 

- m q qq + 

\fl trD.UD^ - l (222) 
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Some explanations about the notation here are in order. Remember that under 
chiral rotations ( Vl,V r ), U transforms like: U —> VrUVl- The unitary matrix U 
is the product of the so-called left and right coset representatives: U = and, 
without lost of generality, one can always choose the gauge where £]) = £# = £. The 
coset representative £, (U = ££,) transforms like: 

£ -> V R £rf = h£Vl h G SU(3)v , (223) 


where h denotes the rotation induced by the chiral transformation (Vl, Vr) in the 
diagonal SU(3)v■ In Eq.( | 222 |) the constituent quark fields Q transform like 


Q hQ, h e SU{ 3) v . 

(224) 

In the presence of external sources^, 


- *(w + - i{vp - O]^} 

(225) 

and 


^ = igDJJt*. 

(226) 


The free parameters of the theory are f n , Mq , and qa- The QCD coupling constant 
is assumed to have entered a regime (below A x ,) where its running is frozen and is 
taken to be constant. 

The merit of this model is that it automatically digests the phenomenological 
successes of the constituent quark model, in a way compatible with chiral symmetry. 
We shall in fact see, that effective Lagrangians of the Georgi-Manohar type, do 
indeed appear in practically all QCD low-energy models where quarks are not 
confined. The weak point of the model is its “vagueness” about the gluonic sector. 
In the absence of a dynamical justification for the “freezing” of the QCD running 
coupling constant, it is very unclear what the “left out” gluonic interactions mean; 
and in fact, in most applications they are simply ignored. 

4.4 Effective Action Approach Models 

The basic idea in this class of models is to make some kind of drastic approxima¬ 
tion to compute the non-anomalous part of the QCD-fermionic determinant in the 
presence of external v [L and a M fields, but with the external s and p fields frozen to 
the quark matrix 

s + ip = M. = diag (m u , rrid, rn s ). 

For this purpose, it is convenient to perform a chiral rotation of the quark fields 
Ql,r{x) = 1(1 ± 75 )q(x) in the initial QCD-Lagrangian: 

Ql{x) =► Q l (x) = f (a;) q L (x), q R =► Qr(x) = ^(x) q R (x), (227) 

§ The original formulation of the model of Georgi and Manohar @ was in fact made without external 
fields. 
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with £ chosen so that ££ = U. Under chiral rotations (Vl, Vr), the quark fields of the 
rotated basis transform like the constituent chiral quarks of the Georgi-Manohar 
Lagrangian: 

Ql,r —> h(x)Q L ,R, (228) 

with h[cp(x)] the rotation in SU( 3)v induced by the chiral transformation (Vl, Vr); h 
is the same object which appears in Eq.( |224|) . The QCD-Lagrangian in the rotated 
field basis, and in Euclidean space has then the following form: 

4cd = ~\ G W G $ + Q d eQ , (229) 

with Dr the Euclidean Dirac operator ( 7 ^ = — are Hermitian Dirac matrices 
with positive metric) 

De = 7 m(<^ + WsG/x + r M — - 75 ^) — ~(E — 75 A), (230) 

where and £ ;/ are the same as in Eqs. (|225|) and ( |22G| ) : and 

s = + £M£, A - CMC - £A4£. (231) 

The E and A terms break explicitly the chiral symmetry. 

The Euclidean effective action Fr(U, v, a, M) associated to Cqc D is then given 
by 

expT E (U, v,a,M) = ^ j T>G^ exp[- J d A z ^ G^G^\ exp T E (G )X ; U , v, a , M) 

= ^ J[VG ft ]expT E (G fi ;U,v,a,M), (232) 

with 


exp T e(G fj,; U, v, a, M.) = J VQVQ exp( J d A z QDrQ ) = det Dr- (233) 

In fact, for the non-anomalous part of the effective action, it is sufficient to consider 
the modulus of the determinant. Then, 

V E (G ll ;U,v,a,M) = ^ log det D ] E D E - (234) 

It is in trying to evaluate the modulus of the fermionic determinant that one 
encounters the first problem. There is need of a regularization; and if possible, a 
regularization which shows the explicit cut-off dependence. The simplest regular¬ 
ization in that respect is the proper-time cut-off method, where 

r E = -~ —Trexp (~tD ] e D e ), (235) 

2 J i/a 2 r 
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where Tr stands for trace over Dirac 7 -matrices, SU (3)-colour matrices, Gell- 
Mann’s flavour-5' U{ 3) matrices, and Euclidean space. 

Since D^ e De is a second-order elliptic operator, the proper-time integrand in 


( |235|) has a well defined power series expansion in powers of r; this is the so called 
heat-kernel expansion: 


[x 


d -tD^„D e 


y) = 


-(x-y ) 2 


I 67 r 2 r 2 


J2H n (x,y)T r 


(236) 


The functions H n (x,y) are known in the literature as Seeley-DeWitt coefficients. 
Here we encounter the second problem : for n > 2, the proper-time integrals are in¬ 
frared divergent. We need an infrared regulator to integrate the large-r behaviour. 
Now, it just happens that, with neglect of the gluonic coupling in the Dirac operator 
in ( p30|) , and for some of the O (p 4 ) terms in the effective action, the proper time in¬ 


tegral in (|235|) is convergent both in the ultraviolet and infrared domains. Therefore, 
the values of the corresponding Lj-coupling constants in this brutal approximation 
turn out to be finite quantities. The results in question are 0: 


SL\ — 4L 2 — Lg — 


N r . 


L‘.\ — T10 — 


N r 


(237) 


48t r 2 ’ “ 1U 96 t r 2 ' 

Numerically, these results, when compared to the phenomenological determinations 
listed in Table 1, are surprisingly good. Any attempt however, to improve on these 
results, and/or to compute other couplings of the low-energy effective Lagrangian, 
necessarily brings in the question of the infrared behaviour of the underlying theory. 

A simple suggestion which has been proposed E3, in connection with the infrared 
behaviour, is to parametrize the effect of spontaneous symmetry breaking, by adding 
to the QCD-Lagrangian a phenomenological order parameter like-term: 

ATqcd = —MQ^RUqL + QlU^Qr), (238) 

which introduces at the same time the [/-held in a way non-invariant under U —> 
—U, and the mass parameter Mq which provides the infrared regulator needed in the 
evaluation of the low energy effective action. It is easy to see that in the presence 
of this term, and withMg > 0 , there is quark condensation S: < > 7 ^ 0 and 

negative. Furthermore, in the presence of this term, it is also possible to evaluate the 
effect of large-iV c gluonic interactions which appear as inverse powers of Mq times 
the appropriate vacuum expectation values of gauge invariant gluonic operators. 
For example, in the chiral limit, and including the leading gluonic contributions, 
one gets ml: 


f 2 = 

J 7 r 


N 

——zAMq 
16t r 2 Q 


log 


A 2 

m 2 q 


+ 


7T 


2 < ^GG > 


6 N r . 


Mq 


+ 


< g 3 GGG > 


360 N r 


m 6 q 


+ ■■■ 


(239) 


There also appear gluonic corrections to some of the previous results for the L~ 
couplings: 


L-a — L 10 — 


N r 


967r 2 


1 + 


7T 


< ^GG > 


5 N r . 






(240) 
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The positive sign of this correction helps towards the agreement with the experi¬ 
mental values. Notice however that the gluon condensate which appears here, is 
an average over configurations from 1/A 2 to roughly 1 /Mq distances; and therefore 
cannot be related easily to the usual gluon condensate which appears in the phe¬ 
nomenology of QCD-sum rules. The predictions for the other L, coupling constants 
to O(jjt). remain the same as those in (|237|). 

Q 

Applications of this approach to the non-leptonic weak interactions have also 
been made. The problem here is to find the effective action of a given four-quark 
operator. For example, we have seen in Sec. 3, that the AS = 2 transitions, after 
the heavy degrees of freedom have been integrated out, are governed by the four- 
quark operator 


Qas =2 = ( s L 'y t "d L )(s L 'y fl d L ), 


modulated by products of flavour-mixing matrix elements, times Wilson coefficient 
functions resulting from the short-distance integration. The Euclidean effective 
action of Qas =2 in an external gluonic background, is given by the expression (see 


Ref. 0 for technical details): 



(241) 


The first term in the r.li.s. is the one induced by the configuration where the two 
currents in Qas =2 are factorized, the term leading in the 1/Ay expansion: the second 
term contains the non-factorizable contributions, which are next-to-leading in the 
1/Ay expansion. These contributions have been evaluated in Ref. 0 including, in 
addition to the well known 0(N 2 ) factorizable contributions, the subleading 0(N C ) 
and 0(a s N c ) terms. The corresponding result for the Bk parameter defined in it 
subsec. 3.3 reads as follows: 



(242) 


The corresponding results for the other coupling constants of the 0(p 2 ) non-leptonic 
effective Lagrangian discussed in Sect. 3, evaluated at the same approximation fJ 
as the B k parameter above, are the following: 



(243) 


and 
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(244) 


+C + (p 2 )- 
+C4(/r 2 ) — 



Ge(p 2 )lbL 5 


N c < ^GG > 
~2 16vr 2 / w 4 


+ 0(a s N c ) 


< > 

~7T~ 


+ 0(a s N c ). 


There are a number of interesting features which emerge from these results, 
worth commenting upon. 


• The results of the so called vacuum saturation approximation (VSA), often used 
in the literature, are obtained from those above when the terms 0(asN c ) 
are dropped. This model calculation shows however that, numerically, the 
neglected VSA-terms are as important as those retained. 


• The results in eqs.( |242D and ( p243| ) satisfy indeed the chiral limit symmetry 
relation between AS = 2 and A I = 3/2, AS = 1 transitions first observed by 
Donoghue, Golowich and Holstein E-l 


Equations ( |242|) , ( 243 ), and ( 244 ) show an interesting correlation between 
the l/AVmorrections to the different couplings. The same correction which 
decreases the G' + - modulated terms, and hence the B[< parameter as well, 
increases the C- modulated term; and hence the A I = 1/2-transitions. It 
can be shown, in full generality, that this is in fact a general property of the 
full 0(A r c )-corrections in the chiral limit El. 


The term proportional to 6/ (jr) in Eq.( |244|) is the result of the next-to- 
leading 1/N C calculation we have already discussed in the subsubsec. 3.2.2. 


• The cancellation of the /^-dependence of the other Wilson coefficients with the 
bosonization of the corresponding four-quark operators, is more involved. To 
exhibit this cancellation explicitly, requires the knowledge of the full dynamics 
of next-to-leading order in the 1/Ay expansion. It is clear that, in as far as 
one doesn’t know the origin of the phenomenological term in eq.( |238|) , this 
cannot be shown. It is possible however to show that the logarithmic ay- 
corrections to the bosonized operators are weighted by the same anomalous 
dimension factors as those of the original four-quark operators; i.e., 


< Q, >~ - log(^j) <Qj>. 


Further calculations within the framework of the Effective Action Approach Model 
discussed here can be found in references @S'" ld 0. A possible generalization of 
the constituent mass ansatz term in ( p38|) to a non-local form has also been sug¬ 
gested 0. Phenomenological applications using a non-local constituent mass term 
can be found in Ref. 0. 
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4.5 The Extended Nambu and Jona-Lasinio Model (ENJL— 
model) 

Since the early work of Nambu and Jona-Lasinio il, there have been many sugges¬ 
tions in the literature proposing models of the type first discussed by these authors, 
as relevant models for low-energy hadron dynamics. [For a recent review where 
earlier references can be found see 0.] The scenario suggested in Refs. US, which 
I shall follow here, assumes that at intermediate energies below or of the order of 
the spontaneous chiral symmetry breaking scale A x , the leading operators of higher 
dimension which, after integration of the high frequency modes of the quark and 
gluon fields down to the scale A x , become relevant in the QCD-Lagrangian, are 
those which can be cast in the form of four-fermion operators, i.e., 


t-'QCD => £qcd + £s,p + £v,a + • • •, (245) 

where 

£s,p = E(^)(4^), (246) 

x ^id 

and 

£v,a = --^5-Gy + L <r+ R\. (247) 

X 'f'lj 

Here i.j denote u. d. and s flavour indices and summation over colour degrees 
of freedom within each bracket is understood; q^y = |(1 ± 75 )q. The couplings 
G sy are dimensionless functions of the ultraviolet integration cut-off A. They are 
expected to grow as A approaches the critical value A x , where spontaneous chiral 
symmetry breaking occurs. (This is the reason why the operators Cs,p and £\/.a 
become relevant.) In QCD, and with the factor N^ 1 pulled out, both couplings G 5 
and Gy are 0(1) in the large- Ay limit. These constants are in principle calculable 
functions of the ratio A/Aqcd- I n practice however, the calculation requires non- 
perturbative knowledge of QCD in the region where A ~ A x , and we shall take G 5 
and Gy, as well as A x , as independent unknown parameters. The ;y index in Cq CD 
means that only the low-frequency modes A < A x of the quark and gluon fields are 
to be considered from now onwards. 

Notice that in QCD, couplings of the type C$,p an< 4 Ty,n a ppear naturally from 
gluon exchange between two QCD colour currents. Using Fierz rearrangement, one 
has in the large- N c limit: 


g 2 sT,a(qY^q)(q i^q) 


N c AT 


CX-s^c 


Eij(?R9Li)(gi,9ffi) 


1 87r 2 a s N c 

N c 7 r 




R]i 


i.e; Gy = G^/4 = j n this case. The two operators Cs.p an d £-v,a have however 
different anomalous dimensions, and it is therefore not surprising that Gg 7 ^ 4Gy 
for the corresponding physical values. 
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If furthermore, one assumes that the relevant gluonic effects for low-energy 
physics are those already absorbed in the new couplings G s and Gy, then 


£qcd => n Pq 

in Eq.( |245| ) with J/) the Dirac operator given in Eq. ( |219|) , where now the gluon field 
G v plays the role of an external colour field source. There is no gluonic kinetic term 
any longer. 

As is well known from the early work of Nambu and Jona-Lasinio 0, the operator 
Cs,p , f° r values of G s > 1, is at the origin of the spontaneous chiral symmetry 
breaking. This can best be seen following the standard procedure of introducing 
auxiliary held variables to convert the four-fermion coupling operators into bilinear 
quark operators. For this purpose, one introduces a 3 x 3 auxiliary held matrix M(x) 
in flavour space; the so called collective held variables, which under chiral ,5'f/(3) 
transform as 

M V R MVl 

and uses the functional integral identity: 

exp [i (<fx -^^-G s J2(VR ( lLj)(qL(lRi)} = 

J JV c ly x Lj 


JVM exp [i Jd A x{ — (q L M\ R + h.c.) — N c 
By polar decomposition 


A 2 1 

(248) 


M = 


with ££ = U unitary and H hermitian. 

Next, we look for translational-invariant solutions, which minimize the effective 
action; 


<9r 


eff 


dM 


H=<H>=Mq ,£ > =l;v=a=s=p=0. 


o. 


The minimum is reached when all the eigenvalues of < H > are equal, i.e., < H >= 
Mq 1; and the minimum condition leads to 


Tr(x 



—2 MqN c 


A x 1 

87T 2 G^ 



(249) 


The trace in the l.h.s. of this equation is formally proportional to < i/ji/j >. The 
calculation however requires a regularization, with A, the ultraviolet cut-off. We 
choose the proper time regularization. [See e.g., Ref. for technical details.] Then 

N M 2 

<i u >= -j^ 4M Q r( .~ L - a 2 * ; (250) 
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and the minimum condition in Eq. ( |249| ) leads to the so-called gap equation: 
Mq 


M 2 q. 


G< 


= M q { exp (——) 


A 2 

x 


Ml , Ml' 

°-r(o, -3-) 


A 2 

x 


The functions 


—f)-/ 


r°° dz 




e~ z z n ~ 2 : 


A 2 

x 


71 = 1, 2, 3,, 


(251) 


are incomplete gamma functions. Equations ( |2 5 0|) and (|251| ) show the existence 
of two phases with regards to chiral symmetry. The unbroken phase corresponds 
to the trivial solution Mq = 0, which implies < TpTp >= 0. The broken phase 
corresponds to the possibility that the coupling G s increases as we decrease the 
ultraviolet cut-off A down to A x , allowing for solutions to Eq. ( |251| ) with Mq > 0 
and therefore < A A > 7 ^ 0 and negative. In this phase the Hermitian auxiliary field 
H(x) develops a non-vanishing vacuum expectation value, which is at the origin of 
a constituent chiral quark mass term [see the r.h.s. of Eq. (|248|j ]: 


—Mq^lU^ qn + qiiUqL ) — —MqQQ , 


like the one which appears in the Georgi-Manohar model; and like the one proposed 
in the effective action approach of Ref. il. 

In the presence of the operator Cv,A- we need two more auxiliary 3x3 complex 
field matrices L^ix) and R^ix) to rearrange the Lagrangian in ( |245|) into an equiv¬ 
alent Lagrangian which is only quadratic in the quark fields. Under chiral (Vl- Vr) 
transformations these collective field variables are chosen to transform as follows: 


L tl V L L^VI , R,, V R R^Vl . 

Then, the following functional identity follows: 

expH / dA x Y%~ Gv ^ ^L^QL^iqi^QLi) + L^ R]}) = 

c 1 X i,j 

JVL.VR, exp[i Jd*x{q L YL,q L + N c ^-^-MrLM, + L ~ R}} . (252) 

It is convenient to trade the auxiliary field matrices L /( (x) and by new vector 

field matrices 

W^ = £L^±?R^, 

which transform homogeneously under chiral transformations (V^, Vr); i.e, 

W& hW^h\ 

H ' h 1 ' 
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with h the SU( 3)y rotation induced by (Vl,Vr). The fermionic determinant can 
then be obtained using standard techniques, like for example the heat kernel ex¬ 
pansion we described earlier. When computing the resulting effective action, there 
appears a mixing term between the fields IF/ 1 and One needs a new redefinition 
of the auxiliary field VF/ ': 


K~ ] - FF + (1 - Sa)H„ , 

in order to diagonalize the quadratic form in the variables and It is this 

mixing which is at the origin of an effective axial coupling of the constituent quarks 
with the Nambu-Goldstone modes: 

a term like the axial coupling which appears in the Georgi-Manohar model, but 
with a specific form for the axial coupling constant g J 4 : 


9 a 


1 


4 A/f2 

1 + Gy^(r(0, 


Mq \ 

A 2 > 
iV x 


(253) 


In terms of Feynman diagrams this result can be understood as an infinite sum of 
constituent quark bubbles, with a coupling at the end to the pion field. These are 
the diagrams generated by the Gy four-fermion coupling to leading order in the 
1/Ay expansion. The quark propagators in these diagrams are constituent quark 
propagators, solution of the Schwinger-Dyson which is at the origin of the gap 
equation in ( |251| ). In the limit where Gy = 0, gA — 1; but in general @ gy 7 ^ 1 to 
leading order in the 1/Ay expansion. 

Kinetic terms for the auxiliary field variables are also generated by the functional 
integral over the quark fields Q and Q. The resulting Lagrangian, after wave- 
function rescaling of the auxiliary fields, has the form of a constituent chiral quark 
model, with scalar S(x), vector V(x). and axial-vector A{x) field couplings: 


y™" = iQY(3„ + r„ - -fir-vM - m q qq 

V aJv 

+ 1;9aQi51^{£» “ - ^-QS(x)Q 

2 JA As 

+ ^tr[d^Sd»S - M 2 s SS] 

- ^tr[(d^V u - d v V li )(&*V v - ff'V^) - 2 MyV^] 

- \tr[{d,A v - dvA^A” - d"A») - 2 M\A^} 

+ + C>(p 4 )terms, (254) 
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where r /; and are the same as those defined in Eqs.( |225|) and ( |226| ), and the 
coupling constants and masses are now expressed in terms of only three input pa¬ 
rameters. As input parameters, we can either fix: G s, Gy, and A x ; or the more 
physical parameters: 

Mq, Ay, g A . 

The coupling constants are then: 


N, 


Si = 0 , Mi/Al) 


(255) 


and the masses: 


and 


/J = I^5 r(0 ’ M eA). 


f\= N ° 2 ~ 2 


1 6 7t 2 3 [3r(0, M q/^x) ~ 2r(1 ’ ^Q /A x } ]’ 


y2 = N c 2 


Ml = 6 Mq—^ , 

~9a 

M a = 6M Qi_g A r(i,M^/A|)> 

r(o,A/g/A|) 


M ^~ 4M Q x 2 r(M^M|) 


3 r(0,A/^/Aj) 


(256) 

(257) 

(258) 

(259) 

(260) 

(261) 


In the absence of the vector and axial-vector four-fermion like coupling i.e., when 
Gy = 0: (]a — 1, My —> oc and Ma —> oo. Then the vector and axial-vector 
interactions decouple, and the model becomes equivalent to the Constitutuent Chiral 
Quark Model of Georgi and Manohar, with g A = 1 and a non-trivial coupling to a 
scalar field. 

The functional integration over the quark fields and the auxiliary S(x), V(x). 
and A(x) fields results in an effective action among the Nambu-Goldstone boson 
particles, with all the couplings fixed by the three parameters Mq, A x , and g A - The 
explicit results one gets for the L t constants which appear in the large - Ay limit at 
0(p 4 ) in the chiral expansion are shown in Table 2. The reason why the constant 
L 7 does not appear in this Table is that, phenomenologically, this constant gets a 
large contribution from the integration of the heavy singlet rf particle. However, 
in the chiral limit, the mass of the rj' is induced by the axial -7/(1) anomaly, which 
only appears to next-to-leading order in the l/A r c -expansion. By definition, the 
ENJL-model, as formulated here, ignores this effect. In order to take these next- 
to-leading effects in 1/7V C systematically, together with the chiral expansion, one 
has to resort to a 1/(3) x 1/(3) formulation of the effective theory 0. The constants 
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Li and Lq are of next-to-leading order in the l/iV c -expansion; this is the reason 
why they do not appear in Table 2 either. We also show in Table 2 the numerical 
results of the fit 1 discussed in Ref. These results correspond to the set of input 
parameter values : 


M q = 265 MeV, A x = 1165 MeV, g A = 0.61. (262) 

Table 2 : The Lj-coupling constants in the ENJL-model of Ref.il, with g A 
defined in Eq.( |253|) , and T n = T(n, Mq/A 2 X ). The second column gives the results 
corresponding to the input parameter values in ( p62|) . The third column gives the 
experimental values of Table 1. 


The Li couplings of 0(p 4 ) in the ENJL-model 

Fit 1 

Experiment 

U = - g\)T , + 2 g\T 2 ] 

0.85 

0.7 ±0.5 

L 2 = 2 L\ 

1.7 

1.2 ±0.4 

L'i - 1671-2 8( [(1 9a) To + 4^(1 g A )^i+ 

-Istpr, - 4r 2 + 3y(r„ - r,) 2 ]} 

-4.2 

-3.6 ± 1.3 

- U] 

1.6 

1.4 ±0.5 

Ls = i£^i[r 0 - fr\] 

0.8 

0.9 ±0.3 

U = Sp\i(l-g 2 A)To + 2g 2 A T 1 ) 

7.1 

6.9 ±0.7 

£io =- Jx)ro + fflO] 

-5.9 

-5.5 ±0.7 


The overall picture which emerges from this simple model is quite remarkable. 
The main improvement with respect to the results obtained in the effective action 
approach model discussed in subsec. 4-4 is 011 the constants L 5 and L 8 , where the 
combined effect of the vector and scalar degrees of freedom leads to rather simple 
results modulated by powers of the g A constant, which agree very well with the 
phenomenological determinations. One of the characteristic features of the ENJL- 
model, is that it interpolates successfully between pure VMD-type predictions and 
those of the constituent chiral quark model. A nice illustration is the result for L 9 
in Table 2 , where the first term is the one coming from vector-exchange, while the 
second one comes from the chiral quark loop integral. 
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There is no difficulty to reproduce the anomalous Wess-Zmnino-Witten func¬ 
tional within the ENJL-model nO 

QCD two-point functions, beyond the low-energy expansion, have also been 
evaluated in the ENJL-model £3. This involves calculations to leading order in the 
l/iV c -expansion (i.e., an infinite number of chains of fermion bubbles; but no loops 
of chains,) and to all orders in powers of momenta Q 2 / A 2 . As a result, vector and 
axial-vector correlation functions have a VMD-like form, but with slowly varying 
couplings and masses. For the transverse invariant functions for example, the results 
are 


4 1} (q 2 ) 


2f 2 mM 2 m 

M 2 m-Q* 


(263) 


and 


n (i)^ = vim , vimMim 

A W ’ Q 2 M 2 (Q 2 ) - Q 2 


(264) 


where 


N r 1 

fv(Q 2 ) = J q dxx(l - *)r(0, X Q = [M Q 2 + x{l - x)Q 2 ]/ Aj]). 
The product 


2 ti(Q 2 ) M v(Q 2 ) = ^ c —— 


' 8 tt 2 Gy 


is scale invariant. 
With 


9a{Q 2 ) = 


l + Gy^/ 0 ]dxV(0,x Q y 


the other couplings are fixed by 


rim=gimtim, 


and the relations 


fymM 2 m = fimMim+fim , ( 265 ) 

and 

fv(Q 2 ) M v(Q 2 ) = fl(Q 2 )M 4 A (Q 2 ). (266) 

The last two equations are the Q 2 -dependent version of the 1 st - and 2 nd -Weinberg 
sum rules 

In the case of the scalar two-point function there appears a pole in the () 2 - 
summed expression at Ms = 2Mq. The case of the other two-point functions is 
somewhat more involved because they mix through the four-fermion interaction 
terms. The corresponding results can be found in Ref.0. Calculations of the low- 
energy behaviour of some three-point functions in the ENJL-model have also been 
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made Corrections due to possible four-quark operators of higher dimension 
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involving derivative terms, have also been studied recently 

In principle, the ENJL-model can be applied to obtain a systematic calcula¬ 
tion of the low-energy constants of the weak non-leptonic Lagrangian, like the 
.Bif-parameter, g 8 and g 2 7 - The difficulties are mainly technical; but progress is 
underway. 


5 Rare Kaon Decays and CP—Violation 

This last lecture is dedicated to a discussion of the theoretical status, within the 
Standard Model, of the CP-violation parameters e and e' which we have introduced 
in Sec. 1. and also to the study of the decay Kl —> n°e + e~ . Other than the 
K —f 7T7T decays which we have discussed, this rare decay mode seems to be the 
most promising candidate to observe direct CP-violation in K physics in the near 
future. 


5.1 The Parameters e and e' Revisited 


In Sec. 1 . we have seen that CP-violation in K —»■ tttt decays is governed by the 
parameters e and e', [cf. Eqs.([52|) and ([001.) ■ | We have also seen that, to a good 
approximation, these parameters can be written as follows: 


e ~ J_ e W4 , Im M 

y/2 V Am Red 0 / ’ 


(267) 


and 


e ~ 


v/2 


3 i(i52-(5o+7r/2) 


Redo /Imdo Imdr 


Redn V Redo Redr 


(268) 


The question we want to discuss here, is the present status in the determination of 
the various components which appear in these expressions. 


5.1.1 The mass difference Am. 

The Kl — K$ mass difference Am is well determined experimentally i: 

Am = m KL - m Ks = (3.522 ± 0.016) x 10“ 12 MeR. (269) 

As we have shown in Sec. 1 ., Am is related to the real part of the off-diagonal 
K° — K° mass matrix, [cf. Eq.(|75j)1: 

Am ~ 2ReMi 2 . 

In the Standard Model, there are both short-distance and long-distance contri¬ 
butions to ReMi 2 - The short-distance contributions, from the effective AS = 2 
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Hamiltonian described in subsec. 3.3 , are modulated by the same unknown four- 
quark matrix element as the contribution to ImM 12 ; i.e., the so called B K paramet er 
dehned in Eq.( |2 13| ) . If one was able to calculate the long-distance component to 
ReMi 2 , it would be possible to fix the value of the B[< factor from the measured 
value of Am. Unfortunately, as I shall next explain, this is not the case. 

In principle, the long-distance contribution to ReMi 2 , and hence to Am, can be 
evaluated in yPT. It appears at the one-loop level, from diagrams with two vertex 
insertions of the effective AS = 1 chiral Lagrangian of 0(p 2 ) described in subsec. 3.2, 
in the presence of the strong interactions. These chiral loops are however divergent, 
and therefore they necessarily bring in new (Dip 4 ), AS = 2 local couplings which 
are not fixed by chiral symmetry requirements alone. Here again, one has to resort 
to models of the low-energy effective action to make further progress. Estimates of 
the long-distance contributions to Am, show that they are important 0. This is 
not surprising: the fact that nominally they are of higher order, Ofp 4 ) in the chiral 
expansion, is here largely compensated by the relative A/ = | enhancement of the 
two lowest 0{p 2 ) vertices. In the absence of a reliable way to calculate the needed 
couplings of the local 0(p 4 ) AS = 2 effective Lagrangian, there appears no way to 
get rid of the unknown Refactor in the ratio : and this is the reason why, for 
phenomenological purposes, one is forced to use the experimental value of Am in 
Eq.fl267|). 

5.1.2 The phase-shift difference <5 2 — 5 0 . 

Contrary to the situation concerning Am, the re — re phase shift difference <f 2 — 5q 
is poorly known experimentally. (It is not easy to extract re — re elastic scattering 
amplitudes from physical observables.) Here, yPT does better! One can extract 
this phase shift difference from the calculation of the elastic re — re amplitudes to 
O(p 4 ), without introducing new unknown parameters, with the result ELJ: 

5/= 2 °(M|) - = -45° ± 6°. (270) 


5.1.3 ReA 0 and ReA 2 . 

There is very little we can do as yet, theoretically, concerning Re/lo and Re/1 2 . 
As we have seen in Sec. 3., their calculation to lowest order in yPT, involves the 
couplings g 8 and g 2 7 , which are not fixed by symmetry requirements alone. The 
models discussed in Sec. 4. are not yet developed to the required degree of precision 
for theory to be useful here. With neglect of electromagnetic corrections, the ratio 
EsiL, can be extracted from the physical branching ratio ’ as a l r eady 

discussed in Sec. 1., [cf. Eqs. (|57|) to (|7C|),1 with the result 


ReA 0 

ReA 2 


+22.2. 


(271) 
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Since Im/1/ <C Reff/, we can also extract Re/lo from the experimental Y(K —► mr) 
rates with the result 

|ReR 0 | = 3.3 x 10 _7 GeV. (272) 

To make predictions about e and e' we still need to know ImM 12 and Im/1 / for 

I = 0 , 2 . 


5.1.4 The e-parameter. 

It is presently known that, experimentally^ : e' <§C e. On the other hand, the Bell- 
Steinberger inequality we have discussed in subsec. 1.2, when restricted to the 2 tt 
intermediate states, can be written as follows: 

( FS ^ Fl + *Am) 2Ree = e r 5 + e'[r 5 (±-) - 2r s (00)]. (273) 

The second term in the r.h.s. is suppressed in two ways: i) because e' <C e; ii) because 
of the A I = 1/2-rule. Furthermore, using the empirical facts that F/ -C F.s' and 
Am ~ r^/2. we arrive at the constraint 


(1 + i)Re e ~ e . 


(274) 


With e given by Eq.(|78|) and e by Eq.( |267| ), this constraint implies the neglect of 


ImAp 

Re/lp 


versus Iu ^ 2 . We conclude that, to a good approximation, e is completely governed 
by the size of M 12 , evaluated in the usual phase convention of flavour mixing in the 
Standard Model. With versus In ^ 12 neglected, it is reasonable to neglect as 
well the long distance effects contributing to ImM 12 . In fact, it has been argued 
that the two effects largely cancel each other To this approximation, e is then 


entirely given by the local structure of the AS = 2 box diagrams which generate 
the four-quark operator in ( |207| ). 

It has become conventional in the phenomenology of flavour dynamics, to use an 
approximate form of the flavour mixing matrix; the Cabibbo-Kobayashi-Maskawa 
matrix^ : 

/ 1-4- A A A 3 ae~ tS 


V = 


-A 

V AA 3 (1 — ae lS ) 


A! 

2 

-AX 2 


1- A 


AX 2 

1 


0(A 4 ), 


(275) 


where A denotes the Cabibbo angle, 


A ~ |Ks| = 0.2205 ±0.0018; 


(276) 


and A and cr are parameters of order one. The phase 6 is at the origin of the CP- 
violation in the Standard Model. The merit of this parametrization is that it clearly 


I See the lectures of Sunil Somalwar EA 

II See the lectures of Roberto Peccei 12 
update of the Unitarity Triangle see Ref. 


for a description of the relevant experiments, 
lot-further details on this parametrization. For a recent 
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shows the orders of magnitude of the various mixing matrix elements. Often, the 
notation 

ae ~ iS = p - irj , (277) 

is also used. In this parametrization, the Jarlskog invariant which governs all the 
observables which violate CP-invariance 0, is given by the product: 

| J CP | = 2 r/X 6 A 2 . (278) 


The A- parameter is fixed from 5-decays into states containing charm: 

A ~ A“ 2 |Kb| = 0.82 ±0.12; (279) 


and a is fixed, once 14b is measured: 

, I Kb 

Act ~ 


V r 


= 0.08 ±0.02 


(280) 


cb 


In terms of this parametrization of the mixing matrix, and to the approximation 
that we have finally adopted, the e-parameter can then be written as follows: 


e = 


—‘° W4 — 1 


G \m 2 w x B k x 


y/T Am 3 Jk k 2M k 2tt 2 ' 
r]X 6 A 2 [r]sS(r c , r t ) - ViS(r c ) + A 2 A 4 (l - p)r/ 2 S(r t )]. 

In this expression, B/< denotes the scale invariant B K parameter: 

Bk = cv s (p 2 )~ 2 / 9 B K (p 2 ); 


(281) 


(282) 


and S(r c ,r t ) and S(r t ) are the box-diagram functions of the ratios r„ = -S-, with 

M w 

rriq the corresponding quark mass running at the pole value: 


S(r t ) = 


n 


1 + 


9 


6 


and 


6r 2 log r t 
1 - r t (1 - r t ) 2 (1 - r t ) 3 

3 r t 


S{r c ,r t ) = r c [log(—) - 1 (1 ± 

r c 4(1 - r t ) 




l°g r t)}- 


(283) 


(284) 


l - n 

The rji parameters correct for short-distance gluon exchanges in the AS = 2 box- 
diagrams. c£a With B factored out, they are renormalization scale invariant, with 
the values** : 


r/i = 1.10, r/ 2 = 0.57, r/ 3 = 0.36. 


(285) 


**The 771-parameter is rat her sensitive to m c and A qcd, which explains the different values found in 
the literature. Here we use !£!! A±± = ( 240 ± 90 )MeP, which corresponds to A±L = ( 350 ± 150 )MeV] 
for the charm quark pole mass we take 0 m c = ( 1.47 ± 0 . 05 )A/eP. 
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Equation ( |281|) , with the value of e fixed from experiment [| e |= (2.26±0.02) x 10 3 ], 
and for a given value of Bk, determines then a hyperbola in the (p,rj) plane. A 
useful approximation to Eq.( |281|) , for quick numerical estimates is 


e 


(3.4 x 10 - 3 )7]A 2 B k 


1 ± 1.3A 2 (1 - p) 



(286) 


5.1.5 Calculations of the B K -parameter. 

There exist several calculations of this parameter using very different techniques. 
The problem is that, although the errors from some of the calculations are rather 
small, the central values are still too dispersed. Two simple calculations are the ones 
we have already discussed in subsec. 3.3: the vacuum saturation approximation LA, 
which gives 

(Bk)vs = 1; (287) 

and the large-A/ prediction, which corresponds to 

B k \n c ->oo = 3/4 . (288) 


None of these two calculations can distinguish between Bk and Bx(p 2 )- 

We also pointed out in subsec. 4-4 that, as first observed by Donoghue, Golowich 
and Holstein Eil, there is a symmetry relation in the chiral limit between AS = 2 
and A I = 3/2, AS = 1 transitions. Using the known K + —> 7r + 7r° decay rate, it is 
then possible to fix Bk in that limit, with the result 


B 


K 


xPTto 0(p 2 ) 


= 0.37. 


(289) 


The calculation of Bk within the effective action approach model we have discussed 
in subsec. 4-4 shows how to reconcile this result with the large-A/, result above. The 
next-to-leading corrections in the 1/Ay-expansion appear to be large and negative. 
The same pattern appears in the 1/Ay approach developped by Bardeen, Buras and 
Gerard. From their matching of short-distances to long-distances, these authors 
give the estimate US: 

B k = 0.70 ± 0.10 . (290) 

This result includes the effect of the one loop 0{p A ) chiral corrections as well, which 
are large and positive. Higher 0(p 4 ) chiral corrections, including the effect of the 
local 0 (p 4 ) terms, evaluated within the effective action approach model, have also 
been recently estimated 0, with the result 


B k = 0.42 ± 0.06 


(291) 


There are a variety of QCD Sum Rules that have been used to estimate Bk as 
well. The description of the techniques involved in these calculations goes beyond 
the scope of these lectures. I shall however give some results. The most recent 
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estimate of Bk using two-point function sum rules of the AS 
operator E3 gives 


B k 


2p.f. sum rules 


0.39 ±0.10. 


2 four-quark local 
(292) 


The same approach reproduces rather well the observed K + —» 7r + 7r° decay rate. 
(It overestimates by less than 15% the coupling constant g 2 7 ■) 

The calculations based on QCD Sum Rules for three-point. functions give a large 
variety of outputs. I shall only quote the result of Ref. L!±l wherefrom one can trace 


earlier results: 


B 


K 


3p.f. sum rules 


0.5 ±0.1 ±0.2 


(293) 


Lattice-QCD simulations of the Bk -parameter tend to find results in the higher 
range of the estimates we have presented here. They are discussed by Steve Sharpe 
in his lectures at this TASI school. 

In my opinion, the only safe claim that theorists can make at present is that the 
value of Bk is likely to lie within the range: 


0.35 < B k < 0.80 . 


(294) 


5.1.6 The ratio e' /e. 


As discussed in the lectures of Sunil Somalwar l__l , there is experimental information 


on this ratio of parameters from the measurement of the ratio of branching ratios 
[cf. Eqs.®,© and ©]: 


with the results: 

fte(-) 

e 


]+_ 

Voo 


|~ 1 ± 6Re(- 


(23 ± 6.5) x 1(T 4 CERN - NA31, 
(7.4 ± 6.0) x 10 -4 FERMI lab - E731. 


(295) 


(296) 


The phases in the expressions of e and e' being practically the same, the theoretical 
prediction for the ratio — turns out to be a real number. If furthermore, we use 
the experimental input in Eq.( |271|) as well as the experimental determination of 
|e| ~ 2.3 x 10~ 3 , we get 


{_ = _1_1_1_ Im ^o (1 _ 22 2 Im/1 ± 

e y/2 (j^ ~ 22.2) (|e| — 2.3 x 10~ 3 ) ReA 0 ’ ImA 0 


(297) 


The dominant contribution to ImA 0 originates in the diagrams which give rise to 
the Penguin-operator Q e . They bring the factors V u \, V c aV* s , and hence the CP- 
violation phase which contributes to Iin/1 0 . We have seen in subsec. 3.2, that to 
lowest order in yPT, the amplitude A 0 takes the form, 


4> = -2| KdlC (& + g±’) V2f n (M, 
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The operator Q e transforms like an (8^,1^) operator and contributes only to gg. In 
the large- N c expansion, its contribution to g§ is next-to-leading. As we discussed in 
3.2.2, it is quite remarkable that, at this order of approximation, this contribution 
can be calculated exactly in terms of known parameters!. Remember the result [cf. 

Eq. (12031)1: 


Imgg 


—ImC6(/r 2 ) 16L.5 


< > 

— ff~ 


{1 + 0(1/N C ) + 0(p 2 )}. 


(298) 


With this result in hand, we can rewrite Eq.(297) in the following way: 


e 


x 


1 1 r)\ 4 A 2 

71 2T2 (2.3 x IQ" 3 ) 


(0.09 ±0.01) 

< V’V’ > 

(Re[g 8 + gJ /2) ] - 5.l) 

p 

L J 7T J 


/ ~ 22 ' 2 Im)p) l 1 + ° (1W + °<±1 ' 


(299) 


Several comments concerning this result are in order: 


• The factor (0.09 ± 0.01) in the r.h.s. is the value of the calculated short- 
distance Wilson coefficient line's, once the overall flavour factor //A 4 A 2 has 
been pulled out. The error here comes from the uncertainties in the top-quark 
mass and in A^r. 


As we already mentioned in subsec. 3.2, the largest uncertainty from the cal¬ 
culated matrix element of the “penguin” Q 6 -operator is due to the poorly 
known value of the quark condensate < ijjip >• [Cf. Eq.( [205|) .l It is cus¬ 
tomary to trade this factor by the value of the strange-quark mass, using 
the Gell-Mann-Oakes-Renner relation in Eq.( |105| ) of lowest order yPT. I pre¬ 
fer not to do that, because then one should also discuss corrections to the 
Gell-Mann-Oakes-Renner relation, which complicates things even further. 


• Because of the large coefficient -22.2- modulating ImA^/ImAo in the r.h.s. of 
the final expression for e'/e above, a reliable estimate of the ratio hii/G/Im/1 0 
is clearly needed. There are two sources of contributions to this term: 

i) Isospin breaking effects, 

ii) Electroweak “penguin” effects, where the gluon exchange of the ordinary 
“penguin”-like diagram is replaced by a photon or a Z vector' boson. Be¬ 
cause of the isovector component in the coupling of the photon and the Z 
to the quark-currents, the electroweak “penguin” diagrams lead to effective 
operators which can induce A/ = 3/2-transitions. 

Both effects i) and ii) have been shown to go in the same direction, and they 
decrease the leading effect of the Q e operator that we have calculated above. 
The effect ii) has become particularly important because, as the top-quark 
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108 


114 


for 


mass rn t increases, for m t > Mw , it- grows roughly as mf. [See Refs, 
a detailed discussion.] Isospin breaking effects due to 7r° — rj — rf mixing have 
been estimated Eds to contibute 


/ lmA 2 
\Iim4 0 


~ 1.4 x 10“ 2 . 

v+v' 


(300) 


• Some attempts to estimate the other O(jj-) and the other 0{jp 2 ) corrections 
in the r.h.s. of Eq.( [299| ) have also been made. The results however are con¬ 
troversial. New efforts in that direction are urgently needed, if we want to 
compare usefully with the next generation of e'/e experiments. 

Due to the various uncertainties we have discussed, it is difficult to claim, at 
present, a theoretical prediction for e'/e better than the limits: 

/ 

10 -4 < - < 2 x 10- 3 . (301) 

5.2 The Decay Kp —> n°e + e~~ and CP—Violation 

Transitions like K —> 717 , with 7 a real photon, are forbidden by gauge invari¬ 
ance. They are allowed, however, for virtual photons 7*, which can produce real 
lepton pairs. Because of the absence of strangenes changing neutral currents in 
the Standard Model, transitions like K —> 7 rl + l~, with l = e,p, are then governed 
by the interplay of weak non-leptonic and electromagnetic interactions. To lowest 
order in the electromagnetic coupling constant they are expected to proceed, in 
principle dominantly, via one-photon exchange. This is certainly the case for the 
K ± —■> 7 r ± l + l~ and Kg —■► decays. The transition K 2 —> 7 r° 7 * —> 

via one virtual photon, is however forbidden by CP-invariance. It is then not 
obvious any longer whether the physical decay K Jj —> 7 will be still domi¬ 
nated by the CP-suppressed y*-virtual transition; or whether a transition via two 
virtual photons, which is of a higher order in the electromagnetic coupling but CP- 
allowed, may dominate. With the possibility of reaching branching ratios for the 
mode Kl —> " (> e + e r as small as 10~ 12 in the near future dedicated experiments, the 
theoretical study of this mode has become of major importance. 


5.2.1 The one-photon exchange amplitude 

The interesting thing about this transition amplitude is that it gets a direct contribu¬ 
tion coming from the terms in the AS = 1 short-distance Hamiltonian proportional 
to the operators Q n and Q i 2 in subsec. 3.1. These are the operators induced by 
“penguin”-like diagrams, where the gluon exchange of the ordinary “penguin”-like 
diagram is replaced by a photon or a Z vector boson. Much the same as in the 
case of the Q fi operator which we have already discussed, these diagrams bring in 
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the factors V u \ VJ*. V cx \ 1%, and hence the CP-violation phase. The mo st, r ecent, cal¬ 
culation of the corresponding Wilson coefficients, can be found in Ref. CEI. Another 
interesting issue concerning this transition, is that the hadronic matrix element here 
is the one of a quark bilinear operator; and by isospin symmetry is related to the 
well known hadronic matrix element of charged decays. The main sources of 
errors in the calculation of this direct transition amplitude are then the mass of the 
top quark, and the CP-violation ^-parameter [see Eq. fl281|) ,1 which modulates the 
overall amplitude. Letting the Bk parameter be within the range in Eq. ([294| ); and 
for a top mass 150GeV < m t < 200GeV, the expected branching ratio from this 
direct source of CP-violation is 


1.3 x l(T i2 < Br (K l -> vr u e + e“ 


I7 dir. 


< 5.8 x 10 


-12 


(302) 


Let me remind you that the present experimental upper limit is 11111 (90% C.L.) 


Br (K l —> 7r°e + e ) < 5.5 x 10 9 , 

exp. 


(303) 


still quite far away. 

The other source of CP-violation, usually called “indirect”, is the one induced 
by the K ^-component of the Kl state, which brings in the CP-violation admixture 
parameter e, i.e., the parameter e to a very good approximation. The problem here 
is then reduced to the evaluation of the CP-conserving transition Kg —> 7r°e + e~. 

The analysis of K —> 7 xl + K decays in general, within the framework of ;\PT 
has been made in Refs. HE3. To (Dip 4 ) in the chiral expansion, the corresponding 
decay amplitudes get contributions both from one chiral loop graphs, and from tree 
level contributions from local operators of O (p 4 ). In fact, only two local operators 
of the O {p 4 ) non-leptonic effective Lagrangian contribute to the CP-conserving 
amplitudes: 

£rff =1 (>) = ~ieG 8 y 2 F 1 ' 1 ' {witr(QA 6 _ i7 E M E i ,) + w 2 tr(Q£ At A 6 _ i7 £ l/ )} , (304) 

J 7T 

where C ll: (x) is the 3x3 flavour matrix field we introduced in Eq.( |185| ); Q the electric 
charge matrix: diag(2/3, —1/3, —1/3); F l ‘ v the electromagnetic field strength tensor; 
A 6 _ i7 the SU( 3) Gell-Mann matrix: \ 6 —i\ 7 ; and w 7 2 are two dimensionless coupling 
constants not fixed by chiral symmetry requirements alone. The overall constant 
G 8 denotes the combination of couplings 

G 8 = KdKs g8 , |G 8 | ~ 9 x 10~ 6 GeV- 2 . (305) 

The mode K + —> 7r + e + e _ turns out to be particularly interesting because both 
its branching ratio and the pion energy spectrum have now been measured in the 
BNL-E777 experiment 0. In full generality, one can predict the K + —> 7 t + 1 + F 
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decay rates as a function of the unknown scale invariant combination of coupling 
constants 


W+ = --(47T 2 )[wi - w£ + 3(\V 2 - 4 L T g)] - - log 


M 


I< 


ml 


(306) 


where w] 2 and Ljj are renormalized couplings at the scale v\ LJ, is the same strong 
coupling constant which governs the mean squared radius of the pion [see the dis¬ 
cussion in section 2.2.2], Viceversa, the measured branching ratio IHZ 


Br (K + -> vr + e + e“) = (2.99 ± 0.22) x 10 


-7 


(307) 


determines two possible solutions for the combination of constants w + . The degen¬ 
eracy between the two possible solutions can be lifted from the measurement of the 
invariant mass distribution of the final lepton pair, (the q 2 dependence; or the 7r + 
energy-spectrum in the K 1 rest-system). A fit to the high q 2 spectrum from the 
latest BNL-E777 data kid favours the positive solution, with the result: 


w + 



0.24 
0.14 ’ 


(308) 


while the positive solution, extracted from the measured decay rate in the same 
experiment corresponds to the value: 

wf=lT (309) 


The two numbers are consistent with each other within less than two standard 
deviations; a fact which provides an independent check of the level of accuracy of 
0(p 4 )-yPT to describe this process. 

Unfortunately, the determination of the combination of constants w + is not 
enough to predict the K$ —> 7T°e + e~ decay rate. To the same order in the chiral 
expansion, the corresponding transition amplitude brings in the combination of 
constants: 

1 1 M 2 

w, =--(47T 2 )[wi - w^] --log—f . (310) 

Clearly, one has to resort to models to go any further in making a prediction. 

In Ref. hJ it was pointed out that the relation w 2 —4L 2 = 0 holds for the divergent 
parts of the corresponding regularized coupling constants. It is also valid for the 
full couplings if, as it happens in many models, the mesonic effective strangeness¬ 
changing current which couples to the virtual photon is required to transform as 
a pure S77(3) octet. In general however, this current is allowed to have terms 
which transform as 10 and 10 S'f/(3)-representations as well. Therefore, as empha¬ 
sized in El the constraint: w 2 = 4Lg, is not required by general chiral invariance 
arguments alone. It is nevertheless tempting to see what the prediction for the K$- 
mode is in the class of models which satisfy this relation. Using the center value 
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for w + — 1.2 from the measurement of the A' + -mode branching ratio quoted above, 
one gets 


Bt(A's —> vr°e + e ) 


~ 5.4 x 10 


-10 


(311) 


which in turn, corresponds to an “indirect” CP-violation branching ratio 


Br (K l 


t r°e H 


l7ind. (wj=4Lg) 


~ 1.6 X 10 


-12 


(312) 


on the lower range of the “direct” CP-violation prediction. 

It is important to analyze the sensitivity of this result to models which do not 
satisfy the constraint w) = 4Lg. This has been done in Ref. t-J and more recently 
in Ref. 1113. The outcome is that the K$ branching ratio is rather sensitive to small 
variations of the octet-dominance constraint. As a result, “indirect” CP-violation 
branching ratios comparable, if not bigger, than the “direct” prediction cannot be 
excluded for the time being. Once more, we find the need to develop good models of 
the low-energy effective action, if one wants to make further progress. There is not 
much else that ,\PT can do here, except wait for the experimentalists to measure 
the rate of Ks —■ y x 0 e + e~; not an easy task! 


5.2.2 The two-photon exchange amplitude 

The transition K 2 —> 7 r° 7 * 7 * —> 7 r°e + e~ is CP-allowed. The lowest non-trivial 
order calculation of this process in yPT involves at least two loops. It is possible, 
however, to obtain a lower bound to the 27 -exchange rate from the calculation of 
the contribution to the absorptive part of the amplitude A(K[ } —>■ 7 i°e + e~) coming 
from the 77 -discontinuity. 

The general kinematics decomposition of the subprocess K® —> 7 T °77 involves 
two invariant amplitudes. With 

^ 2 (p)^ 7 rV)+ 7 (<?i)+ 7 fe), p 2 = M\ p ' 2 = M' 2 , q\ = q\ = 0, 

gauge invariance and Lorentz invariance restrict the form of the transition amplitude 
to 


OL, 


M[K 2 (p) -> TT (p'h(qih(q2)} = G g — e^q^e^) [A(y,z) (q%q" - q x ■ q 2 g^) 


+2 B(y, z) (p ■ q x q%p u + p ■ q 2 qiP^ - p^p u qi ■ q 2 ~ p ■ qiP ■ q 2 9^) 


where 


y = 


P-(Qi- 92 ) 


and 


z = 


hi + ®) ! 


(313) 


(314) 


M 2 M 2 

Bose symmetry requires the amplitudes A and B to be symmetric functions of q\ 
and q 2 , i.e. 

A(-y,z) =A(y,z), B(-y,z) = B(y,z). (315) 
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To lowest non-trivial order in yPT, only the amplitude proportional to A con¬ 
tributes; and in fact, to that order, A depends only on 0 i.e., the invariant mass 
squared of the 77 -pair. None of the local terms of 0(p 4 ) in the effective AS = 1 
non-leptonic Lagrangian can contribute to this decay. As a result, the full con¬ 
tribution at lowest non-trivial order in the chiral expansion, comes only from the 
finite chiral one loop amplitude. The predicted ^-distribution has a very char¬ 
acteristic shape; being highly peacked at the higher end of the spectrum. This 
prediction Ell has been subsequently confirmed experimentally by the CERN-NA31 
collaboration £29. However the 0(p 4 ) predicted branching ratio: 


Bt(K l —> 7 r°e + e ) = 6.8 x 10 
turns out to be smaller than the observed rates: 


-7 


Br (K l 


7 r°e + e ) = 


(1.7 ±0.2 ±0.2) x 1(T 6 
( 1.86 ±0.60 ±0.60) x KT 6 


NA31 

E731 


120 


121 


(316) 


(317) 


It has been claimed, however that phenomenologically expected higher order effects 
can explain this discrepancy l 122 H 123 l . 

The term proportional to the A amplitude in Eq. ( [> 1 3| ) has a tensor structure 
which, when inserted in the two-body 77 phase space integral to get the corre¬ 
sponding K® —► 7 r°e + e“ amplitude, requires a flip in the helicity of the electron 
line for the transition to be allowed; a fact which brings in the electron mass m e 
as an overall suppression factor in amplitude. This amplitude structure leads to a 
branching ratio: 


Br(AT 2 ° 


7T°77 


7r°e H 


■) 


±77) 


~ 5 x 10 


-15 


(318) 


too small, by several orders of magnitude to be competitive with the one-photon 
exchange, CP-violating amplitude, contributing to K r — > 7T°e + e . 

The interesting issue here is that, contrary to what happens with the contribu¬ 
tion from the A ( 77 )-amplitude which we have just discussed, the tensor structure 
proportional to the R-amplitude in Eq.( [313|) , when inserted in the two-photon 
phase space integral of the process R 2 —> 7 ^ 7 * 7 * —> 7r°e + e _ has contributions 
which are allowed without requiring an electron helicity flip; i.e., they are not m e - 
suppressed. Although the B ( 77 )- amplitude first appears at 0(p () ) in yPT; the fact 
that it can induce a helicity allowed transition to Jt 2 — 1 ► 7 r° 7 * 7 * —> n 0 e + e ~, makes 
it “the dominant amplitude” for this process. 

There have been several estimates in the literature 


HE 


of the B amplilude in 
Eq.( |313| ), based on vector meson dominance ( VMD ) models. Here, however, one 
has to be careful not to overestimate its size and therefore spoil the observed shape 
of the invariant 77 -mass distribution in the Kl —> n 0 e + e~ decay. Many of the early 
models are in fact now ruled out by the NA31 experiment; but they have been very 
useful to sharpen our views on the role of vector mesons in yPT in general, with 
the results 11 we already discussed in subsec. 4 . 2 . 
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Within the framework of yPT, one can parametrize the local 0(p 6 ) contributions 
to K l —> 7 r °77 induced by vector-exchanges, (including vector-exchanges in direct 
weak transitions,) by an effective vector coupling ay, such that 

5 = -2a v G 8 M 2 K -. (319) 

7T 

In terms of this parametrization, the contribution to the —> 7r°e + e _ amplitude 

is as follows: 

M K 2 -»■ 7 T°e + e“) |B(rr)= iG 8 ^-a v u(k )7 • p— (320) 
which leads to a branching ratio 0 


Br (K 


7T°77 


ttV 


■8(77) 


= 4.4 x 10 




(321) 


The NA31-collaboration 
an invariant 77 -mass: M, 
the following limits: 


12 c 


77 


has performed an analysis of Kl —> 7T°e + e events with 
< 240MeV, using this parametrization, and obtained 


0.32 < ay < 0.19 (90% C.L.) 


(322) 


The rate in Eq.( |321|) , and hence the result in Eq.( |322| ) does not take into account, 
however, the non-polynomial structure of the 5-amplitude in Eq. (|313|) due to chiral 
loops. A recent analysis IH1 which tries to fold the phenomenology of both local and 
non-polynomial effects in the 5-amplitude with the observed data from the NA31- 
experiment, results in a branching ratio for Kl —■► 7T°e + e~ from the 77 -discontinuity: 


Br (K l 


7 r°e H 


27 Abs. 


0.13 x 10~ 12 , a v = 0, 


1 s v in - 12 




(323) 


The value ay = —0.9 is the one which reproduces the observed ^-spectrum in the 
decay rate of Kl —> 7 T° 77 , in the presence of the non-polynamial ansatz for the 
5-amplitude which has been used in Ref. LEHl. The two effects combined also raise 
the predicted K L —■> 7r°77 branching ratio to 1.6 x 10~ 6 , in good agreement with the 
experimental values in Eq.( |317|) . 

As a summary, the theoretical prospects for more refined predictions on the 
decay Kl — > 7r°e + e _ are the following: 


• The branching ratio expected from “direct” CP-violation is rather well known. 
The uncertainty in Eq.( [302|) will be reduced once the top quark mass is better 
determined; and the 5#-factor is pinned down more accurately - either by 
theoretical improved calculations, or phenomenologically- . 
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• The error in the branching ratio expected from the CP-conserving transition 
induced via two intermediate photons, can be reduced with a combined effort 
of more accurate measurements of the mode Kl —> vr °77 on the one hand, 
and an improvement in the phenomenological ansatz of the B amplitude in 
the theoretical analysis on the other. 

• The largest uncertainty, at the moment, comes from the “indirect” CP 
violation transition. Here, the only way I can see to make progress is via 
the development of good models of the QCD low-energy effective action; e.g., 
the extension of the succesful ENJL-model in the strong sector to non-leptonic 
weak transitions. Parallel to this theoretical effort, there should be, of course, 
some progress as well in obtaining more experimental results in rare K decays, 
so as to have enough observables to test the models. 
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